Critical Casimir forces between homogeneous and chemically striped surfaces 
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Recent experiments have measured the critical Casimir force acting on a colloid immersed in a 
binary liquid mixture near its continuous demixing phase transition and exposed to a chemically 
structured substrate. Motivated by these experiments, we study the critical behavior of a system, 
which belongs to the Ising universality class, for the film geometry with one planar wall chemically 
striped, such that there is a laterally alternating adsorption preference for the two species of the 
binary liquid mixture. For the opposite wall we employ alternatively a homogeneous adsorption 
preference or homogeneous Dirichlet boundary conditions. By means of mean-field theory, Monte 
Carlo simulations, and finite-size scaling analysis we determine the critical Casimir force acting on 
the two parallel walls and its corresponding universal scaling function. We show that in the limit of 
stripe widths small compared with the film thickness, on the striped surface the system effectively 
realizes Dirichlet boundary conditions, which generically do not hold for actual fluids. Moreover, 
the critical Casimir force is found to be attractive or repulsive, depending on the width of the stripes 
of the chemically patterned surface and on the boundary condition applied to the opposing surface. 
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I. INTRODUCTION 

As an intriguing consequence of their presence, fluctu- 
ations of an embedding medium may manifest themselves 
in terms of effective forces acting on its confining bound- 
aries. The critical Casimir force is such a fluctuation- 
induced force which arises due to the emergence of long- 
ranged thermal fluctuations if a fluid close to a second- 
order phase transition is confined between surfaces. This 
phenomenon, first predicted by Fisher and de Gennes [l| 
is the analogue of the Casimir effect in quantum elec- 
trodynamics 0. Reference Q provides a recent review 
which illustrates analogies as well as differences between 
these two effects and guides the reader towards further 
reviews of the subject and the pertinent original litera- 
ture. 

The dependence of the critical Casimir force on the 
distance between the confinements and on temperature 
is characterized by a universal scaling function, which is 
determined by the bulk and surface universality classes 
(UC) 4, 5] of the confined system. It is independent of 
microscopic details of the system, and its form depends 
only on a few global and general properties, such as the 
spatial dimension d, the number of components of the 
order parameter, the shape of the confinement, and the 
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type of boundary conditions (BC) 

In recent years the critical Casimir effect has attracted 
numerous experimental (9l4l7j and even more theoreti- 
cal investigations. Critical Casimir forces can be inferred 
indirectly by stud ying wetting films of fluids close to a 
critical end point [I8lfl9|. In this context, 4 He wetting 
films close to the onset of superfluidity Q and wetting 
films of classical [10l ] and quantum [ill ] binary liquid mix- 
tures have been studied experimentally. Only recently 
direct measurements of the critical Casimir force have 
been reported [T3 - [l6l | by monitoring individual colloidal 
particles immersed into a binary liquid mixture close to 
its critical demixing point and exposed to a planar wall. 

Not only the strength of critical Casimir forces can be 
tuned by small temperature changes but even their sign 
depends on the BC of the confining boundaries. The two 
interfaces of a 4 He film impose a symmetry-preserving 
Dirichlet BC (denoted by (o)) on the superfluid order- 
parameter at both sides of the film, which causes attrac- 
tive critical Casimir forces leading to a thinning of the 
film near the A-transition [§, HH, 03 ■ However, for clas- 
sical binary liquid mixtures (or simple fluids), surfaces 
preferentially adsorb one of the two species of the mix- 
ture (or the gaseous or the liquid phase of a simple fluid, 
respectively). This corresponds to symmetry-breaking 
BC (denoted as (+) or (— ) BC) acting on the order pa- 
rameter which is, e.g., the concentration difference in a 
binary liquid. Within the theoretical description (±) BC 
are realized by surface fields and the (o) BC by their 
absence. 
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The emergence of long-ranged thermal fluctuations 
close to a second-order phase transition leads to a meso- 
scopic extent of the adsorption layer close to surfaces 
with (±) BC. Depending on whether the adsorption pref- 
erences of the confining surfaces of the fluid are the 
same (±,±) or different (+,—), critical Casimir forces 
acting on them are either attractive (±, ±) or repulsive 
(+,-) [S OS El- The critical Casimir force be- 
tween walls with (±) BC is the combined effect of the 
change of the fluctuation spectrum due to the confine- 
ment and the interference of the adsorption layers, which 
are present even within mean-field theory. The shapes of 
the adsorption layers themselves are strongly influenced 
by non-Gaussian fluctuations, i.e., they differ from mean- 
field predictions. In this sense, the effective forces act- 
ing on surfaces which confine a (near-) critical fluid pro- 
vide a classical analogue of the Casimir effect both in the 
case of symmetry-breaking and in the case of symmetry- 
preserving BC. 

Early on theoretical investigations of the critical 
Casimir force used, to a large extent, field-theoretical 
methods (see, e.g., Ref. [20| for a list of references). Only 
recently Monte Carlo (MC) simulations have allowed for 
their quantitatively reliable computation. Early numer- 
ical simulations for the critical Casimir force have been 
employed in Ref. |2l| for the film geometry with later- 
ally homogeneous BC. More recently the critical Casimir 
force has been determined by numerical simulations for 
the XY UC 22-28] , which describes the critical proper- 
ties of the superfluid phase transition in 4 He, as well as 
the Ising UC [2(| [H, 0, [H H-fll] which describes, inter 
alia, the experimentally relevant demixing transition in 
a binary liquid mixture. 

Since Casimir forces may affect or empower future de- 
vices on the micro- and nano-scale, their modifications 
due to the presence of nano- or micro-structures on the 
substrates has been a topic of intense research during the 
last decade. Recent theoretical and experimental studies 
of QED Casimir forces (see, e.g., Ref. [35 1 an d references 
therein) as well as critical Casimir forces 36| for topologi- 
cally structured substrates exhibit remarkable deviations 
from the corresponding ones for planar walls as well as 
the occurrence of lateral forces. However, only chemi- 
cally patterned substrates allow for interesting combina- 
tions of attractive and repulsive critical Casimir forces 
so that among the various realizations of the critical 
Casimir effect, the force in the presence of a chemically 
patterned substrate has recently attracted particular in- 
terest [H,|37j]. 

Experiments with binary liquid mixtures as solvents 
have been used to study critical Casimir forces acting on 
dissolved colloids close to a chemically structured sub- 
strate |14l - ll6l |. which creates a laterally varying adsorp- 
tion preference for both components of the solvent. Such 
kind of systems have been investigated theoretically for 
the film geometry within mean- field theory [38j . within 
Gaussian approximation [39j . and with MC simulations 
in a three-dimensional film geometry in the presence of a 



single chemical step [20|. The critical Casimir force in the 
presence of a patterned substrate has also been studied 
in the case of a sphere near a planar wall within the Der- 
jaguin approximation [I(|[4lj]. If the lateral chemical pat- 
terns do not consist of stripes with sharp chemical steps 
between areas of strong but opposite adsorption prefer- 
ences, one faces spatial regions characterized by surface 
fields of medium strength. This case has been studied so 
far only for laterally homogeneous BC in the presence of 
variable boundary fields. This case gives already rise to 
interesting crossover phenomena, which have been stud- 
ied within mean- field theory |42l , by exact calculations in 
two spatial dimensions [43l.l44l|. and with MC simulations 

Mm 

Motivated by the aforementioned experimental results, 
and based on previous investigations by two of the au- 
thors [20j], here we present a MC study of a three- 
dimensional lattice model in the film geometry, repre- 
senting the Ising UC in the presence of a chemically 
striped substrate. Moreover, we compare the universal 
scaling functions of the critical Casimir forces obtained 
from these MC results with the corresponding mean-field 
results, which we obtain by generalizing a previous study 
38] and which are valid in d = 4 spatial dimensions. 
We employ periodic boundary conditions in the lateral 
directions and different BC for the two surfaces confin- 
ing the slab. To this end, we consider a film of thickness 
L confined along the normal z-direction on one side by a 
surface at which the order parameter of the fluid exhibits 
a laterally homogeneous BC which corresponds either to 
strong adsorption (+) or to the so-called ordinary surface 
transition (o) [HQ. The other side of the film is confined 
by a surface which is periodically patterned by stripes 
leading to strong, alternating adsorption preferences cor- 
responding to (+) or (— ) BC, respectively, varying along 
the lateral ^-direction. 

Here, we focus on stripes of equal width S+ = 5_ = 
P/2 corresponding to half of the period P along the re- 
direction, so that the important geometrical parameter 
is given by n = S+/L, which relates the width of the 
stripes to the film thickness (see Figs. Q] and [2]). Within 
the lattice model this system is realized by either fixing 
the Ising spins in the upper surface to +1 or imposing an 
open boundary by not fixing them, whereas the lower sur- 
face consists of alternating stripes of equal width, where 
the spins are fixed to +1 and —1. The chemical steps 
separating the stripes are taken to be sharp. 

Our results show that in the limit of stripe widths small 
compared to the film thickness, the lower surface effec- 
tively realizes Dirichlet BC. Thus the system reduces for 
k — > to (+, o) or (o, o) BC, and, in order to be able to 
compare with this limiting case, here we also consider a 
film in which both surfaces have a laterally homogeneous 
BC from the outset (see Figs.[3]and|4|). This may provide 
a novel possibility of studying also symmetry-preserving 
BC for simple fluids and binary liquid mixtures which 
are difficult to establish experimentally otherwise [l6l ]. 

In order to extract universal quantities from MC simu- 
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FIG. 1. Film geometry confined by a laterally homogeneous 
upper surface and by a lower surface with alternating stripes 
of equal width. At both surfaces the spins are fixed. We 
choose S+ — S- so that the period P = S+ + S- = 2S+. 




FIG. 2. Film geometry confined by an upper surface with 
open BC and by a lower surface with alternating stripes of 
equal width with fixed spins. 

lations, it is important to take corrections to scaling into 
account in order to be able to extrapolate data for sys- 
tems of finite size L to the thermodynamic limit L — > oc. 
In particular, in the standard three-dimensional Ising 
model, scalin g co rrections are proportional to L - ^, with 
cj = 0.832(6) [451 ]. The presence of non-periodic bound- 
ary conditions, such as in the direction normal to the 
film, gives rise to additional scaling corrections, the lead- 
ing one being proportional to L , which is numerically 
difficult to disent ang le from the previous one. Following 
Refs. [H, [H [H Hi H2, in order to avoid the simulta- 
neous presence of these competing corrections, we have 
studied a so-called improved model [46], for which the 
leading scaling corrections oc L~ u are suppressed for all 
observables so that the correction oc L _1 becomes the 
leading one. 

This paper is organized such that in Sec. [IT] the finite- 
size scaling behavior, as expected for the system under 
study, is established. In Sec. IHII we introduce the lattice 
model studied here. In Sees . II VI and IVl we present our MC 
results for the critical Casimir force at T = T c and for 
the universal scaling function of the critical Casimir force 
at T 7^ T c , respectively. The corresponding results ob- 
tained within mean-held theory (d = 4) are presented in 




FIG. 3. Film geometry confined by a laterally homogeneous 
upper surface with fixed spins and by a lower surface with 
open BC. 
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FIG. 4. Film geometry confined by a lower and an upper 
surface both with open BC. 



Sec. I VI I and compared with the actual behavior in d = 3 
in Sec. lVIll We summarize our main findings in Sec. lVIlIl 
In App.[X]we provide certain important technical details 
of the MC simulations. In App. [B] we report details of 
the determination of the bulk free energy density which 
is needed in order to compute the critical Casimir force. 



II. FINITE-SIZE SCALING AND CRITICAL 
CASIMIR FORCE 

In this section we recall the finite-size scaling (FSS) 
behavior of a system in the film geometry L x lA~ 1 in 
d spatial dimensions, which in the thermodynamic limit 
exhibits a second-order phase transition at the temper- 
ature T = T c . Here, we restrict ourselves to the BC 
described above; a broader discussion of finite-size scal- 
ing for non-periodic BC can be found in Ref. [2(|. In the 
following, for the sake of brevity, we do not analyze sep- 
arately the FSS behavior of the BC illustrated in Figs. [3] 
and |H where there are no stripes. These two cases can 
be obtained by taking the limit k = S+/L — > in the 
BC of Figs. [1] and [21 respectively. 

In the critical region and in the absence of an exter- 
nal bulk field, the free energy density T per fc^T of the 
system (i.e., the free energy divided by Lin -1 fe^T) can 
be decomposed into a singular contribution and a non- 
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singular background term: 

T{t, L,L b S+)= -F (s) (t, L,L b S+)+ -F (ns) (t, L,L {{ ,S+), 

where t — (T — T c )/T c is the reduced temperature. The 
non-singular background .F( ns ) can be further decom- 
posed into specific geometric contributions, correspond- 
ing to bulk, surface, and line contributions. According to 
renormalization-group (RG) theory 47] and neglecting 
corrections to scaling, in spatial dimension d the singular 
part of the free energy density obeys the following scaling 
property: 

^(t,L,L h S+) = ^f(T,K,p), 



k = S+/L, 
P = L/L h 



(2) 



where v is the critical exponent of the bulk correlation 
length and £q~ is its non-universal amplitude: 



(3) 



The function f[r,n,p) is a universal scaling function, 
i.e., it depends only on the bulk universality class and 
on the BC applied at the two surfaces. As in Ref. [11], 
the scaling ansatz in Eq. ^ generalizes the one for lat- 
erally homogeneous BC by an additional dependence on 
the scaling variable k. In the following we neglect the de- 
pendence on the aspect ratio p = L/Lu because here we 
are interested in the film geometry with Lu L. In this 
limit and for the BC considered here, the dependence on 
the aspect ratio is expected to be negligible. Our MC 
data support this observation (see also the discussion in 
Sec. [TV] below). The bulk free energy density /bulk(*) is 
defined as 



/ bu lk(t)= lim T(t,L,L h S+) (4) 

L,L\\— too 



and it is independent of the BC. Analogously to Eq. (TTJ, 
fbuik(t) can also be decomposed into a singular contribu- 
tion and a non-singular background: 

A»ik(t) = &(*) + (5) 

with /^ k (* -> 0) cx \t\ du = \t\ 2 ~ a , where a is a standard 

bulk critical exponent. The excess free energy fe$ is 
defined as the remainder of the free energy density J-^ 
after subtraction of the bulk contribution: 

/W(t, L, L h S + ) = .7*0 (i, L,L h S+) - /«k(t). (6) 

According to Eq. ([2]) it exhibits the following scaling be- 
havior: 
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f^(t,L,L b S+) = J? A[T = t(L/^) 1 ,k = S + /L 



(7) 



The critical Casimir force Fc per area L 
UbT is defined as 



(d-i) 



Fr 



d(LfW 



dL 



t.Lu ,S.i 



and per 



(8) 



Due to Eqs. the critical Casimir force exhibits 

the following scaling behavior: 
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F C (t, L, L h S+) = ^9(r = t {L/i+) 1/v , k = S+/L 

(?) 

where 0(r, n) is a universal scaling function. At the crit- 
ical point one has r = 0, so that at criticality the force 
is given by 



Fc(t = 0,L,L h S + ) = I3 Q(K), (10) 



with 



G(k) = (9(0, k). 



(11) 



In the limit of very small stripes, i.e., k — > 0, the char- 
acter of a striped surface effectively approaches the one 
for a homogeneous one with (o) BC. Thus, the scaling 
functions of the critical Casimir force approach the ones 
for the critical Casimir force acting on two homogeneous 
surfaces with (+,o) or (o, o) BC, respectively, i.e., 

K ^ (6,+ At), (+) vs. stripes for L > S+, 

V+/o{T,k) > < 

I "(o,o)(t), (o) vs. stripes for L ^> b+, 

where the subscript + jo indicates the corresponding type 
of BC at the homogeneous surface. 

On the other hand for very broad stripes, i.e., n — > oo, 
the limiting behavior for the case of a homogeneous (+) 
wall opposite to a striped surface (Fig. [TJ is given by 
the average of the two homogeneous cases for (+, +) and 
(+, — ) BC, respectively. In this case, i.e., for n 1 
the system effectively corresponds to the one for isolated 
chemical steps opposite to a homogeneous wall, connect- 
ing regions which are almost laterally homogeneous and 
correspond to (+, — ) or (+,+) BC. Thus, for fixed film 
thickness L, the contributions from the nearly isolated 
chemical steps to the critical Casimir force per unit area 
vanish as cx P -1 , so that the corresponding contributions 
to the scaling function vanish as a k" 1 . The asymptotic 
behavior for L <C >5+ of the universal scaling function 
for the critical Casimir force for a (+) wall vs. a striped 
surface is therefore given by 



9+(t,k > 1) 



)(r) 



->M) 



(13) 



where the proportionality constant is related to the ef- 
fects induced by the chemical steps separating the stripes 
and may depend on r. 

Similarly to Eq. (TTBl . for the case of a (o) wall vs. a 
striped surface (Fig. [2]) 0(t, k) approaches 



9 g (t, k > 1) - 0( +)O )(t) oc k \ 



(14) 
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because 0(+, o )(t) = ( _ o) (t). 

For r < 0, due to the presence of the chemical steps 
between the stripes, interfaces form, which separate the 
domains of positive and negative order parameter. As 
will be discussed below, for the case of a (+) wall oppo- 
site to a striped surface as well as for a (o) wall opposite 
to a striped surface and n < 2, these interfaces align on 
average parallel to the film surfaces. By contrast, for a (o) 
wall opposite to a striped surface and k > 2 the emerg- 
ing interfaces for r < preferentially align perpendicular 
to the film surfaces in order to minimize the interface 
area. As discussed in Sec. IVII below, for the latter case 
the proportionality constant in Eq. ((13)1 is determined 
by contributions from these interfaces and is given by 
-R^t^, where R a = a (to ) d ~ 1 / (k B T c ) is the univer- 
sal amplitude ratio for the interfacial tension a = tro I * I M 
associated with the spatially coexisting bulk phases and 
fj, = (d — l)v is its critical exponent. Thus, for the limit 
r <C — 1 and n > 2 the scaling function of the critical 
Casimir force between a (o) wall and a striped surface 
approaches 

6 (t « -1, K > 2) ~ c+ , o) (t) - —\t\^. (15) 

Accordingly, the limits for r — >• — oo and K->oodo not 
commute. 

III. LATTICE MODEL AND OBSERVABLES 

In order to compute the critical Casimir force for a bi- 
nary liquid mixture close to its critical demixing point, 
as in Ref. I20| w e study the so-called improved Blumc- 
Capel model |4§L l49l ] as a representative of the 3D Ising 
universality class. It is defined on a three-dimensional 
simple cubic lattice, with a spin variable S{ on each site 
i which can take the values Si — — 1, 0, 1. The reduced, 
dimensionless Hamiltonian for nearest neighbor interac- 
tions is 

« = -0 SiSj+Dj^Sl Si = -1,0,1, (16) 

<ij> i 

so that the Gibbs weight is exp(— H) and the partition 
function is 

Z08,L,L||)=2exp(-W), (17) 

{C} 

where {C} is the configuration space of the Hamiltonian 
given in Eq. (|16p . We note that the partition function 
in Eq. p7|) depends implicitly also on the BC (see the 
discussion below). In line with the convention used in 
Refs. 0, S Ea, HqI, in the following we shall keep 
D constant, considering it as a part of the integration 
measure over {Si}, while we vary the coupling parame- 
ter /3, which is proportional to the inverse temperature, 
/3 ~ 1/T. In the limit D — > — oo, one recovers the usual 
Ising model, because in this limit any state for which 



there is an iq such that Si — is suppressed relative 
to the states {Si — ±1}. For d > 2, the model exhibits 
a phase transition at j3 c — (3 C (D) which is second order 
for D < Z?tri and first order for D > D tl i. The value of 
D tr i in d = 3 has been determined as D tl i = 2.006(8) in 
Ref. [5l|, as D tT i ~ 2.05 in Rcf. 52], and more recently 
as D t ri = 2.0313(4) in Ref. 

We consider a three-dimensional simple cubic lattice 
L z x L x x L y , with L y = L x and periodic BC in the lateral 
directions x and y. For the two confining surfaces we 
employ the BC shown in Figs. [TJU The BC illustrated in 
Fig. [U are realized by fixing the spins at the two surfaces 
z = and z = L z — 1, so that there are L z — 2 layers 
of fluctuating spins. The spins at the upper surface z = 
L z — 1 are fixed to +1, and the lower surface z = mimics 
a patterned substrate, so that the surface is divided into 
stripes of equal width s+ and alternating BC with the 
spins fixed to +1 or —1, respectively. 

Here and in the following all lengths are measured in 
units of the lattice constant a. The size L z indicates 
the total number of lattice layers, including eventually 
the layers of fixed spins. Therefore the thickness L, the 
lateral size Ln, and stripe width S+ are related to the 
dimensionless lattice lengths L z , L x , and s+ according 
to L = (L z — l)a, L« — L x a, and S+ — s + a, respec- 
tively. For the sake of simplicity, here and in the fol- 
lowings Sees. HVlTVl we employ a slightly different def- 
inition of the scaling variables t and k. We consider 
ti = t(L z /^ l ) 1 / v and k ; ee s+/L z , where §jj = ^/a 
is the dimensionless non-universal amplitude of the cor- 
relation length on the lattice, measured in units of the 
lattice constant. Accordingly, we also redefine the aspect 
ratio as pi = L z / ' L x . By comparing these new definitions 
with the previous ones introduced in Eq. ([2]), we observe 
that, for L -> oo, t(L z /^) 1/v = t{L/^f' v + 0{1/L), 
s+/L z = S+/L + 0(1/ L), and L z /L x = L/L [{ + 0(1/ L). 
Therefore, the FSS limit, i.e., the limit L z — > oo at fixed 
t\, ki, as well as the limit of vanishing aspect ratio pi — > 0, 
are unaltered by these new definitions. In order to avoid 
a clumsy notation, in the following we omit the index /. 

Here we consider the limit of a vanishing aspect ratio 
p = L z / L x — > 0, which is obtained via extrapolation by 
computing the critical Casimir force for three different 
aspect ratios p < 1 (see the discussion in the following 
sections). As discussed at the end of Sec.|Hj for the BC 
illustrated in Fig. [TJ in the limit p — > the subsequent 
limit k = s+ / L z — > oo corresponds to the presence of an 
isolated chemical step. In such a geometry, the isolated 
chemical step gives rise to a line defect which, in turn, re- 
sults into a linear aspect ratio dependence of the critical 
Casimir force. In the limit of vanishing aspect ratio the 
force reduces to the mean value of the force for homoge- 
neous (+, +) and (+, — ) BC, for which the two surfaces 
disp lay the same (resp., opposite) adsorption preference 
[20| (compare with Eq. (fT3"|) ). In the opposite limit k — » 0, 
the lower surface is expected to effectively realize Dirich- 
lct BC (compare the upper part of Eq. (fT2]0 . In order to 
analyze this case, as a reference system we study a film 



6 



geometry L z x L x x L x with periodic BC in the lateral 
directions x and y, fixed spins at the surface z = L z — 1, 
and open BC on the lower surface, so that there are L z — 1 
layers of fluctuating spins. This geometry is illustrated 
in Fig. [3] In the following, we shall denote this BC as 
(+,o). 

In addition we consider the three-dimensional film ge- 
ometry L z x L x x L x with periodic BC in the lateral 
directions x and y, with fixed spins at the lower surface 
z = 0, and open BC at the upper surface, so that there 
are L z — 1 layers of fluctuating spins. For the lower sur- 
face z — we employ a pattern such that the surface is 
divided into alternating stripes of equal width s + with 
the spins fixed to either +1 or — 1. This geometry is il- 
lustrated in Fig. [2] Two interesting limiting cases arise 
from this geometry. In the limit of large stripes, i.e., for 
k = S+ / L z — > oo and for vanishing aspect ratio, the lower 
surface effectively realizes an isolated chemical step. In 
analogy with the results of Ref. [20| , in this limiting case 
the critical Casimir force is the mean value of the force 
for (+,o) and (—,6) BC, which correspond to a film ge- 
ometry where one of the confining surface implements 
Dirichlet BC, and the other surface exhibits a homoge- 
neous adsorption preference for one of two components 
of the fluid. In the absence of an external bulk magnetic 
field these two BC are equivalent. Therefore we conclude 
that in the limit n = s + /L z — > oo and for vanishing as- 
pect ratio, the critical Casimir force for the BC of Fig. [2] 
reduces to the force for the (+, o) BC illustrated in Fig. [3] 
(compare with Eq. (fT4]) ). 

In the opposite limit n — > 0, the lower surface effec- 
tively realizes Dirichlet BC, so that the system reduces 
to a film geometry with Dirichlet BC on both surfaces 
(compare with the lower part of Eq. (|12|l ). In order to 
analyze this limit, as a reference system we consider here 
a three-dimensional film geometry L z with pe- 

riodic BC in the lateral directions x and y and open BC 
at both surfaces, so that there are L z layers of fluctuating 
spins (see Fig. 0}. In the following we shall denote this 
film BC as (o, o). 

For the lattice model corresponding to Eq. (|T6|) . the 
scaling behavior discussed in Eqs. ©, ©, and © is valid 
only up to contributions due to corrections-to-scaling. 
We distinguish two types of scaling corrections: non- 
analytic and analytic ones. The non-analytic corrections 
are due to the presence of irrelevant operators. In this 
case in Eq. © additional scaling field contributions arise, 
which are characterized by negative RG dimensions. In 
the FSS limit, i.e., for L z — > oo, t — > at fixed £/L z , this 
results in the following expression for the singular part 
of the free energy density in the absence of external 
bulk fields: 



F^{t,L = a{L z -1),L,| =aL x ,S+ 



asj 



1 



Lk>l 



(18) 



relevant operators and gi are smooth functions which 
are universal up to a normalization constant. The lead- 
ing correction is given by the operator that has the 
least negative dimension. This is usually denoted by 
u, so that the leading scaling corrections are cc L~ w . 
For the standard three-dimensional Ising model one has 
cj = 0.832(6) [45|. In a family of models characterized by 
an irrelevant parameter A, it can occur that for a certain 
choice of A the amplitude of the leading correction-to- 
scaling term oc vanishes. In these so-called improved 
models, the observed scaling corrections usually decay 
much more rapidly, i.e., as L~ U2 with 0J2 = 1.67(11) 
according to Ref. [54J and uj 2 — 1.8867038 according 
to Ref. [551 ] for the three-dimensional Ising universality 
class. This scenario holds for the Blume-Capel model 
described by Eq. (|16[) . where D is an irrelevant param- 
eter for D < Ari- At D = 0.656(20) (H the model is 
improved. In the present work we fix D = 0.655, which 
is the value of D used in most of the recent simulations 
of the improved Blume-Capel model [H, EH HI El • For 
this value of the reduced coupling D the model is critical 
for /3 = /3 c = 0.387721735(25) g|. 

Another type of scaling corrections is provided by so- 
called analytic scaling corrections, which can stem from 
various sources. Nonlinear terms in the expansion of the 
scaling field r [56[ result in scaling corrections oc L z 1 . 
Analytic corrections can also be due to the boundary con- 
ditions: BC which are not periodic in all directions induce 
additional corrections, which are proportional to L~ , It 
was first proposed in Ref. [57j ] in the context of studying 
surface susceptibilities, that such scaling corrections can 
be absorbed by the substitution L z — > L z + c, where c is 
a non-universal, temperature-independent length. Re- 
cently, this pro perty has been checked numerically in 
Refs. I26l . l58l .l59l] for the XY model with free surfaces, in 
Ref. j29| for the Ising model with homogeneously fixed 
surface spins, and in Refs. [13, Hll for the Ising model 
with laterally inhomogeneous surfaces. 

Here we study the critical Casimir force using the im- 
proved Blume-Capel model according to Eq. (fTB)) . On the 
basis of the above discussion, for such a model the lead- 
ing scaling corrections are expected to be proportional 
to L^ 1 . Furthermore, assuming that also in this case in 
leading order such a scaling correction can be absorbed 
by the substitution L z — » L z + c, Eq. ([9]) is replaced by 



F c (t, L = a(L z - 1), L|| = aL Xl S. 




where y^ < 0, i > 1, are the RG-dimensions of the ir- 



(19) 



In the case of laterally homogeneous BC in Figs. [3] 
and 21 the dimensionless quantity c (such that ca is a 
length) enters only via the volume factor and via the 
scaling variable r. Scaling corrections to Eq. (fT9")l are ex- 
pected to decay as oc L Z U2 (with uj 2 = 1.67(11) [Hi] or 
uj 2 ^ 1.8867038 [H, see above). 

We introduce the reduced energy density 
E(j3, L Zl L x , s + ) in units of — ksT, which is used 
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in order to compute the critical Casimir force: 

E(0,L z ,L x ,s+) = ^{Y1 )> ( 2 °) 



<i3> 



where V = L Z L X is the total number of spins and (. . .) 
denotes the thermal average. (Note that, according to 
Eq. (HU), has no contribution ~ J2i S ?-) The re- 

duced free energy density F(0, L z , L x , s + ) is defined as 



F(0,L z ,L x ,s+) 

1 / Z(/3, L = a(L z - 1), L\\ = aL a 



V lU \ Z(0, L = a{L z - 1), L,| = aL 



(21) 



Thus F(/3, L Zl L Xl s+) is the free energy per spin and 
in units of — fc^T. It is normalized such that F(0 — 
0, L z , L x , s+) = 0. With this normalization one has 

rP 

F(0,L z ,L x ,s+) = d0'E(0',L z ,L x ,s + ). (22) 
Jo 

The relation between F(t, L z , L x , s+) and the reduced 
free energy density F((3,L z ,L x ,s + ) defined in Eq. (121)) 
is given by 

F(0,L z ,L x ,s+) 

= -F(t, L = a(L z - 1), L\\ = ai x , 5+ = as + ) 

+ .F(i — > oo, L = a(L 2 — 1), £ii = aL x , S+ = as + ). 

(23) 

Finally, the reduced bulk free energy density Fbuik(/3) is 
defined by taking the thermodynamic limit of Eq. (|2ip : 



F bulk (/3)= lim F(0, L Z ,L X , s+). (24) 

L Z ,L X — >oo 



IV. CRITICAL CASIMIR AMPLITUDE AT T c 

In order to determine the critical Casimir force at T c , 
we follow the approach introduced in Ref. [24| and also 
used in Refs. [Ill [H, [32[ , which we briefly describe here. 
For two reduced Hamiltonians Hi and H2 associated with 
the same configuration space {C} we construct the con- 
vex combination H(X) 



where (H2 — H\)\ is the thermal average of the observ- 
able H2 — Hi with the statistical weight exp(— H(X)). 
For every A this average is accessible to standard MC 
simulations. Finally the integral appearing in Eq. (|27[) 
is performed numerically, yielding the free energy differ- 
ence between the systems governed by the Hamiltonians 
H2 and Hi, respectively. 

We apply Eq. (|2"?| with Hi as the Hamiltonian of the 
lattice L z x L x x L x with the BC illustrated in Figs. [THU 
and H2 as the Hamiltonian of the lattice (L z — 1) x L x x L x 
plus a completely separated two-dimensional layer of 
non-interacting spins governed by the reduced Hamilto- 
nian of Eq. (TiT)l) with — 0, so that both Hamiltonians 
share the same configuration space. This layer can be 
inserted into the film by varying the coupling (1 — A)/? 
with its neighboring planes between and 0. With this 
we evaluate the following quantity: 



1 f 1 

I (0, L Z ,L X , s+) = -j / d\(H 2 - Hi) 

L x JO 



(28) 



By using the definitions of the excess free energy (Eq. ©) 
and of the critical Casimir force (Eq. ©) one finds [20] 

I(0,L z ,L x ,s + )=F hnlk (0) 



+ Fc ( t, L = a [ L z - - ) , i|| = aL x , S + = as + 



(29) 



where corrections oc L~ 2 have been neglected. In com- 
puting the critical Casimir force, the derivative in Eq. © 
is implemented by a finite difference between the free en- 
ergies of a film of thickness L — a(L z — 1) and of a film of 
thickness L — a = a(L z — 2), so that the resulting critical 
Casimir force corresponds to the intermediate thickness 
a(L z — 3/2). This choice ensures that in the FSS limit 
no additional scaling corrections oc L~ l are generated 
[20l ]. By inserting Eq. into Eq. we obtain the 
following scaling form for I (0,L z ,L x ,Sj r ): 



I(0,L z ,L x ,s+)=F hnlk (0) 

I / (L, 1 




l/v 



L z - 



(30) 



H(X) = (1 - X)Hi + XH2, A e [0, 1] . (25) At the bulk critical temperature Eq. (J3DJ) turns into 



This Hamiltonian H(X) leads to a free energy F(A) in 
units of ksT [601. Its derivative is 



dF(X) _ 2^{c} S>r e 



I (0 C , L Z ,L X , s+) = F bu \ k (0 c ) 

1 V - 



(31) 



dx 



(26) 



Combining Eqs. (|2"5)l and (j2"6")l we can determine the free 
energy difference as 

F(l) - F(0) = J dX^^- = J dX(H 2 - Hi)x, (27) 



(L z -l/2 + cf \L z -l/2 + c 
Equation pip can be rewritten as 

I(0 c ,L z ,L x ,s+) = F hulk (0 c ) 



(32) 



with c ' given by 



2.5 



c 

2 



30 (k) 9k 



(33) 



In a series of MC simulations, we have evaluated the 
quantity / (/3 C , L z , L x , s+) for lattice sizes L z = 8, 12, 
16, 24, 32, and 48 with the BC illustrated in Fig. Q] for 
K = 1/4, 1/2, 1, 2, and 3 as well as with the BC of Fig. H 
which corresponds to the limit k — > 0. We have also com- 
puted / ((3 C , L Z ,L X , s + ) for lattice sizes L = 8, 12, 16, 24, 
and 32 with the BC illustrated in Fig. [2] for k = 1/4, 
1/2, 3/4, 1, 2, and 3 as well as with BC of Fig. H which 
corresponds to the limit n — > 0. Certain important de- 
tails of the simulations are reported in App. [X] Since we 
are interested in the film geometry, which corresponds 
to the limit of a vanishing aspect ratio p = L z /L Xl we 
have simulated every BC for three aspect ratios p < 1/8, 
such that there is always an even number of stripes in 
the lower confining surface. An odd or non-integer num- 
ber of stripes would give rise to a line defect which in 
turn, for p — > 0, would result into an unwelcome linear 
aspect-ratio dependence [20] . Within the present numer- 
ical accuracy, for p < 1/8 the MC data do not show a 
visible dependence on p. Thus we consider our results 
obtained for non- vanishing p < 1/8 as a reliable extrap- 
olation to the limit p — > 0. A posteriori, this also jus- 
tifies the scaling ansatz in Eqs. (fT))- (fl~Tj) . in which the 
dependence on p has been neglected. We have simulated 
the Blume-Capel model with the Hamiltonian given in 
Eq. (fl6|) . choosing the values of the reduced couplings as 
D = 0.655 and (3 C = 0.387721735. This corresponds to 
the critical point of the improved model [45], for which 
the Eq. (f52"|) is expected to describe correctly the correc- 
tions to scaling. We have fitted our MC data directly to 
the quantity / (f3 c = 0.387721735, L z , L x , s+) in Eq. (1521 . 
leaving -Fbuik(/3c), 0, and c ' as free parameters. In order 
to control a possible systematic error due to subleading 
scaling corrections, we have repeated the fits discarding 
the smallest lattices. For the BC of Figs. [1] and O and for 
various values of ratio k, in Tables HI and HTl we report the 
fit results as a function of the smallest lattice size L m i n 
taken into account for the fit. In Tables IIIII and IIVI we 
report the corresponding fit results for the BC of Figs. [5] 
andU 

Inspection of the the fit results tells that we gener- 
ally reach a good x 2 /DOF ratio and the results appear 
to be stable with respect to the choice of L m in- {DOF 
is the number of degrees of freedom, i.e., the number 
of statistically independent points minus the number of 
fit parameters) While there is a clear dependence of the 
Casimir amplitude on k, as expected the critical bulk 
free energy density Fbuik(/3c) does not exhibit a depen- 
dence on k. Furthermore the latter is in agreement with 
the value F hulk ((3 c ) = 0.0757368(4) reported in Ref. [Hj]. 
By conservatively judging the variation of the resulting 
with respect to L m i n , from Tables HI and ITTI we obtain 
the following estimates for the BC shown in Figs. [T] and 

El 
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FIG. 5. Critical Casimir force amplitude 0+(k) = 9+(0, n) 
(see Eqs. © and CCD) at T c for the BC of Figs. [Q and [3] and 
for k = S+/L = 0, 1/4, 1/2, 1, 2, and 3 as inferred from 
Tables Hand |IT] (see Eqs. ^-^). The amplitude at k = 
is obtained for the (+, o) BC illustrated in Fig. [3] The dashed 
line provides a smooth interpolation. This line saturates at 
6+(k -> oo) = (e {+ , +) + e {+ ,_)) /2 = 2.386(5) 20], which 
is indicated by the dotted line. The omitted statistical error 
bars defined as one standard deviation and calculated with 
the standard Jackknife method (see, e.g., Ref. [62]) are com- 
parable with the symbol size. 



(+) vs stripes: 



H 


-(« 


= 0) = 


(+)O ) = 0.492(5), 


(34) 


H 


-(« 


= 1/4) 


= 0.62(1), 


(35) 


H 


-(« 


= 1/2) 


= 0.85(1), 


(36) 


H 


-(« 


= 1) = 


1.383(4), 


(37) 


H 


-(« 


= 2) = 


1.875(6), 


(38) 


H 


-(« 


= 3) = 


2.053(5). 


(39) 



The subscript + indicates the homogeneous (+) BC on 
one of the confining surfaces. These amplitudes are 
shown in Fig. [S] As expected, for decreasing values of k 
the critical Casimir amplitude 0(k) approaches the corre- 
sponding value for (+, o) BC. In particular, 0+(« = 1/4) 
is only 26% larger than + (O). In the opposite limit 
k — > oo, + (k) approaches the critical Casimir ampli- 
tude for a single chemical step: 0+(k — > oo) = 2.386(5) 
;20]. In particular, 0+(k = 3) is only 14% smaller than 
0+(k — > oo). In the whole sampled region, + (k) is a 
positive and monotonically increasing function of n so 
that the critical Casimir force at T c is always repulsive. 
The critical Casimir amplitude 0+(O) = 0(+, o ) f° r (+>°) 
BC can be compared with, e.g., the amplitude 0( + + ) 
resulting from homogeneous BC (+,+), for which the 
two confining surfaces exhibit the same adsorption prefer- 
ence. Within mean-field theory one has Q/+ )0 )/©rj- +) = 
-1/4 [2l|. According to the MC results of Ref. (29|', one 
has ©(+,+) = —0.820(15) so that the ratio between the 
two amplitudes is 0r +jO )/0( +)+ ) = —0.60(1). Thus the 
fluctuations produce a significant dependence of this ra- 
tio on the spatial dimension. Accordingly, one concludes 
that in d = 3 mean-field theory captures only the quali- 
tative behavior of the critical Casimir force. Our result 
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12 



16 



k -» : (+,o) 
X 2 /DOF = 8.7/15 
F bulk (/3 C ) = 0.0757369(2) 
6+ = 0.492(5) 
c ' = 0.36(3) 
X 2 /DOF = 8.0/12 
F hulk (p c ) = 0.0757368(2) 
9+ = 0.495(10) 

c ' = 0.40(9) 
X 2 /DOF = 7.4/9 
F bulk (/3 C ) = 0.0757368(2) 
6+ = 0.50(2) 
c ' = 0.4(2) 



k = 1/4 

X 2 /DOF = 12.3/15 
Fbuik(/? c ) = 0.0757369(1) 
6+ = 0.622(5) 
c ' = -0.48(2) 
X 2 /DOF = 7.5/12 
Fbuik(/? c ) = 0.0757368(2) 
6+ = 0.634(11) 
c ' = -0.39(7) 
X 2 /DOF = 6.5/9 
Fbuik(/3 C ) = 0.0757368(2) 
9+ = 0.63(2) 
c ' = -0.39(15) 



k = 1/2 

X 2 /DOF = 16.1/15 
F bu ik(/3c) = 0.0757375(1) 
6+ = 0.845(5) 
c ' = -0.44(1) 
X 2 /DOF = 13.2/12 
FbuikOSc) = 0.0757375(2) 
6+ = 0.84(1) 
c ' = -0.44(5) 
X 2 /DOF = 7.7/9 
Fbuik(ySc) = 0.0757372(2) 
6+ = 0.88(2) 
c ' = -0.23(12) 



TABLE I. Fit of our MC data at T c for the BC of Figs. H and [T] to Eq. J32) with free parameters F bulk (^ e ), 6+(k = s+/L g ), 
and c '. Lmin is the smallest lattice size taken into account for the fit. DOF denotes degrees of freedom. The quoted error 
bars of the fit parameters correspond to one standard deviation; see, e.g., Ref. [6l[ for a discussion of the method of minimum 
X 2 data fitting. 



12 



16 



K=l 

X 2 /DOF = 8.9/15 
Fbuik(/3 C ) = 0.0757370(1) 
6+ = 1.383(4) 
c ' = -0.264(8) 
X 2 /DOF = 4.8/12 
F bu ik(/3 C ) = 0.0757369(2) 
6+ = 1.387(8) 
c ' = -0.25(2) 
X 2 /DOF = 4.2/9 
F bu ik(/? c ) = 0.0757369(2) 
9+ = 1.394(12) 
c ' = -0.22(5) 



k = 2 

X 2 /DOF = 12.7/15 
Fbuik(/? c ) = 0.0757366(1) 
6+ = 1.875(5) 
c ' = -0.138(8) 
X 2 /DOF = 11.0/12 
Fbuik(/3 C ) = 0.0757367(2) 
6+ = 1.869(8) 
c ' = -0.15(2) 
X 2 /DOF = 7.1/9 
Fbuik(/3 C ) = 0.0757368(2) 
9+ = 1.86(1) 
c ' = -0.18(4) 



k = 3 
X 2 /DOF = 9.0/15 
F bu ik(/3c) = 0.0757369(1) 
6+ = 2.053(4) 
c ' = -0.097(5) 
X 2 /DOF = 7.0/12 
FbuikOSc) = 0.0757369(2) 
6+ = 2.048(8) 
c ' = -0.11(2) 
X 2 /DOF = 5.0/9 
F b uik(y8c) = 0.0757369(2) 
6+ = 2.05(1) 
c ' = -0.09(3) 



TABLE II. Same as Table [Qfor « = s+/L z = 1, 2, 3 and for the BC of Fig. 



12 



16 



« — )■ : (o, o) 
X 2 /DOF = 6.9/12 
F bu ik(/3 C ) = 0.07573678(9) 
6 = -0.030(2) 

c ' = 0.8(2) 
X 2 /DOF = 3.7/9 
F bu ik(/3 C ) = 0.0757368(1) 
6 = -0.030(5) 

c ' = 0.7(7) 
X 2 /DOF = 3.2/6 
F bu ik(/? c ) = 0.0757368(3) 
e o = -0.035(15) 
c ' = 1.5 ±2.3 



k = 1/4 

X 2 /DOF = 7.5/12 
F bulk (/3 C ) = 0.0757369(1) 
9 = -0.039(2) 

c ' = 0.02(9) 
X 2 /DOF = 3.8/9 
F bu ik(/3 C ) = 0.0757370(2) 
e o = -0.045(5) 

c ' = 0.5(4) 
X 2 /DOF = 2.5/6 
FbuitCSc) = 0.0757369(3) 
9 = -0.038(10) 
c ' = -0.2 ± 1.2 



k = 1/2 

X 2 /DOF = 13.7/12 
F bulk (/3c) = 0.0757369(1) 
9 = -0.054(1) 
c ' = 0.07(6) 
X 2 /DOF = 11.0/9 
fWk^c) = 0.0757369(2) 
9 = -0.053(3) 

c ' = 0.0(3) 
X 2 /DOF = 7.5/6 
F bulk (/?c) = 0.0757368(3) 

e = -0.05(1) 

c ; = -0.1(9) 



TABLE III. Same as Table U for the BC of Figs, fj and 
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-fcmin K — 3/4 

8 X 2 /DOF = 9.9/12 

Fbuik(/3c) = 0.07573679(9) 
e o = -0.062(2) 
c ' = 0.37(6) 
12 X 2 /DOF = 7.4/9 

Fbuik(/3 C ) = 0.0757367(1) 
e o = -0.058(4) 
c' = 0.1(2) 
16 X 2 /DOF = 4.4/6 

F bulk (/3 C ) = 0.0757369(3) 
9 = -0.07(1) 
c ' = 0.9(8) 



X 2 /DOF = 8.0/12 
F bulk (/3 C ) = 0.0757370(1) 
e o = -0.032(2) 

c ' = 1.3(2) 
X 2 /DOF = 7.9/9 
Fbuik(/? c ) = 0.0757370(2) 
9 = -0.032(5) 

c ' = 1.2(7) 
X 2 /DOF = 3.4/6 
F bulk (/3 C ) = 0.0757367(2) 
6 = -0.021(8) 
c ' = -0.8 ± 1.7 



X 2 /DOF = 13.1/12 
Fbuik(/3 C ) = 0.0757365(2) 
6 = 0.185(4) 
c ' = 0.34(5) 
X 2 /DOF = 8.9/9 
J Fbuik(/3 c ) = 0.0757369(3) 
e o = 0.173(9) 
c ' = 0.04(20) 
X 2 /DOF = 5.6/6 
F bu i k (/3 C ) = 0.0757361(6) 
O = 0.20(3) 
c ' = 0.9(6) 



X 2 /DOF = 12.0/12 
Fbuik(/3 C ) = 0.0757368(2) 
e o = 0.287(4) 
c ' = 0.36(4) 
X 2 /DOF = 8.3/9 
J Fbuik(/3 c ) = 0.0757366(3) 
9 = 0.292(10) 
c ' = 0.45(14) 
X 2 /DOF = 6.6/6 
-Fbuik(/3 C ) = 0.0757363(6) 
9 = 0.30(2) 
c ' = 0.65(35) 



TABLE IV. Same as Table EU for « = s+/L z = 3/4, 1, 2, 3 and for the BC of Fig. [2] 



for 9+(k = 0) = 6 
the result 9 



(+,o) 



(+.o) = 0.492(5) is in agreement with 
while it is not 



e (+i0) = 0.33 



0.497(3) of Ref. 1311 
compatible with the earlier results [2ll 
and 0.416 obtained with the e— expansion method and 
0.375(14) obtained by MC simulations [2l| . 

Inspecting the results reported in Tables Mil and IIV1 
we obtain the following estimates for the BC shown in 
Figs. M and H 

(o) vs stripes: 



e (n 


= 0) = 


e (o , o) = -0.030(5) 


(40) 




= 1/4) 


= -0.039(6), 


(41) 


e ( K 


= 1/2) 


= -0.053(3), 


(42) 


&o(K 


= 3/4) 


= -0.062(4), 


(43) 




= 1) = 


-0.032(3), 


(44) 


&o(K 


= 2) = 


0.18(1), 


(45) 




= 3) = 


0.287(5), 


(46) 



where the subscript o indicates the homogeneous Dirich- 
let BC on one of the two confining surfaces. These am- 
plitudes are shown in Fig. [5] As expected, for decreas- 
ing values of n the critical Casimir amplitude o (k) ap- 
proaches the corresponding value 6( 0)0 ) for (o, 6) BC, 
while in the opposite limit n — > oo it approaches slowly 
the value 0(+, o ) f° r BC. Moreover, the critical 

Casimir amplitude changes sign: it is attractive for n = 
and repulsive for k — > oo. Inspecting Fig. [6l we can es- 
timate that o (k) vanishes for k ~ 1.2. Remarkably, 
different than 0+(k) in Fig. [SJ the critical Casimir am- 
plitude o (k) is not monotonic, but exhibits a mini- 
mum close at K fa 3/4. Our result for Q (k = 0) = 



9 



(o,o) 



-0.030(5) is in agreement with the recent MC 



result 6 (k = 0) 



9 



(o,o) 



-0.028(16) of Ref. 



and also with the earlier results [21( O o (0) = —0.0278 
and —0.0328 obtained with the e— expansion method and 
9 o (0) = -0.023(4) obtained by MC simulations [H|- 

Finally, we can test the validity of Eq. (|3"3"|) by studying 
the behavior of the scaling corrections in the limit n — » 0. 
To this end, we consider the BC of Fig. [T] and we take 
the limit of n — > at fixed L z , i.e., s+ — > in Eq. (|3"TT) . 
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FIG. 6. Critical Casimir force amplitude Q (k) = o (Q, k) 
(see Eqs. © and dTTJ) at T c for the BC of Figs. [5] and g] and 
for n = S+/L = 0, 1/4, 1/2, 3/4, 1, 2, 3, as inferred from 
Tables ITTT1 and ITVl (Eqs. (gDJ-®). The amplitude at k = 
is obtained for the (o, o) BC illustrated in Fig. [4] The dashed 
line provides a smooth interpolation. This line saturates at 



e 



(+.o) 



0.492(5) (Eq. p^l. which is indicated 



by the dotted line. The comparison with the thin full line tells 
that 6„(k) changes sign at k ~ 1.2. The omitted statistical 
error bars are comparable with the symbol size. 



Assuming that O(k) is analytic close to n — 0, we obtain 



/ Q9 C , L z , L x , s+ ->■ 0) = F bulk (/3 C ) + 



e+(o) 



(L z -l/2 + c)3' 
(47) 

A comparison of Eq. (|47p with Eq. (I3^t gives c '(k — > 
0) = c, a result which could also be obtained by taking 
the limit k —¥ in Eq. f|33|) . On the other hand, in 
the limit s+ — > 0, the system effectively realizes the BC 
shown in Fig. [31 but still in the presence of only L z — 2 
fluctuating layers of spins (as for the BC in Fig. [T] with 
s + > 0). According to the convention fixed in Sec. IIII1 
this corresponds to (+, o) BC for a film with L z — 1 layers 
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and thickness a{L z — 2): 

I (0c L z , L x , s + -t 0) = 7 (+)0 ) (J3 C , L z - 1, La,) 

e (+ , o) (48) 



F bulk (/3 C ) + 



(^-1-1/2 + ^)3' 



where the subscript (+, o) denotes explicitly the BC of 
Fig. [3] with the convention of Sec Mil and where we have 
used Eq. ([3"2"|). By comparing Eq. (gTJ) with Eq. (gS]) we 
finally obtain: 



lim c '(k) 

K— >0 



(+,o) 



1. 



(49) 



We can extract c '(+. Q ) = 0.36(4) from the fit results of 
Table|T]for the (+, o) BC. This result is in marginal agree- 
ment with the result c '(+, ) = 0.42(2) of Ref. [3l[ in 
which the same improved Blume-Capel Hamiltonian as 
the present one has been simulated [63j. Using Eq. (|49|) 
we obtain c = c '(k — > 0) = —0.64(4). Inspecting the fit 
results of Tables HI and [TT1 we observe that c '(k) varies 
smoothly with k and indeed approaches the value of 
c = —0.64(4) for k — > 0. According to the results of 
Eqs. f3"3 ]) -([59" |) and due to Fig.[5J the coefficient multiply- 
ing (c — 1/2) in Eq. (|33p is positive. This would imply 
that, due to c — 1/2 < 0, c '(k) < c. However, within the 
current numerical precision such an inequality appears 
to be not satisfied by the fit results reported in Tables U 
and mi This suggests that the ansatz of Eq. (TK)1) does not 
completely capture the scaling corrections for the striped 
BC. One may need to modify in addition the second scal- 
ing argument of 9 in Eq. (|19p . for example by replacing 
L with L + aN, with N an integer number depending 
on the convention used to measure the film thickness, 
or more generally, by introducing a second non-universal 
length. A similar analysis of the scaling corrections for 
the BC shown in Fig. [5] is beyond the presently available 
numerical precision. 



V. THE CRITICAL CASIMIR FORCE SCALING 
FUNCTION 

The determination of the critical Casimir force off criti- 
cality has been performed using essentially the al gor ithm 
introduced in Ref. [23| and also used in Refs. |26l - [29l |31| . 
By using the definition of the critical Casimir force given 
in Eq. ([8]) , the definition of the reduced free energy den- 
sity given in Eq. (|2ip , and the definition of the reduced 
bulk free energy density given in Eq. (IM1) . the critical 
Casimir force can be expressed as 



Fc \ t, L = a { L z , L« = aL r , S+ = asj . 

V V 2/ " J (50) 

= AF(/3,L z ,L x ,s+)-F hulk (l3), 

where 

AF(/3, L z ,L x ,s+) = L z F(p, L z ,L x ,s+) 

-(L z -l)F(/3,L z -l,L x ,s+). 

(51) 



Analogous to Eq. (|21)1) . in Eq. the derivative in 

Eq. @ is implemented by a finite difference between the 
free energies of a film of thickness L = a(L z — 1) and 
of a film of thickness L — a = a(L z — 2), so that the 
resulting critical Casimir force corresponds to the inter- 
mediate thickness a(L z — 3/2). This choice ensures that 
in the FSS limit no additional scaling corrections oc L z l 
are generated [20j . The reduced temperature t is given 
by t = (& - /3)/ p, with p c = 0.387721735(25) H|. As 
in Eq. (|2"9"]l . in Eq. (|S"0"j) corrections oc L z 2 have been 
neglected. We note that AF(/3, L z , L x , s + ) — > Fb u ik(/3) 
for L Z ,L X — > oo, which is in accordance with the van- 
ishing of the critical Casimir force in the limit of large 
volume. Another useful relation follows from a compari- 
son of Eqs. I|50] l and ([29" ]) : 



(52) 



AF(/3, L z ,L x ,s+) = I OS, L Z ,L X , s+) . 



Instead of using the coupling parameter approach as in 
Sec. IIVI here we compute the free energy differences by 
sampling the internal energy density E(/3, L z , L Xl s+) for 
various values of j3 and for film thicknesses a(L z — 1) 
and a(L z — 2). Then AF(f3, L Z ,L X , s+) is computed by a 
numerical integration of Eq. ((22)) . For doing so, it is very 
useful to observe that it is not necessary to perform the 
integral in full between f3' = and /3' = j3 [29(. In fact, 
by inserting a lower cut-off (3q m to the integral appearing 
in Eq. ([22| one can effectively compute the difference 
between the critical Casimir force and the force at the 
inverse temperature /3q. This implies that the critical 
Casimir force can be expressed as 



F c \ t,L = a yL z - - I , L|| = aL x , S+ = as + 
= AF(L z ,L x ,s+;p,/3 ) - (F bulk (/3) - F bulk (/3 )) 



+F C ( t ,L = a ( L z - - ) , i|| = aL x , S+ = as^ 



(53) 



with 



AF(L Z ,L X , s + ;(3, (3 ) = L z / dp'E(f3' ,L z ,L x ,s+) 
-(L z -1) / dp'E(P',L z -l,L x ,s+), 

Jf3o 

(54) 

and to = (/3c — M/Pa as the reduced temperature cor- 
responding to the lower cutoff /3q. Since for L — a(L z — 
1) 3> C the critical Casimir force vanishes oc exp(— L/£), 
one can neglect the last term in Eq. (|55)l if the corre- 
lation length £ at the lower cutoff j3o is much smaller 
than L = a(L z — 1). Moreover, due to Eqs. (|52l) and 
((50"]) the last term in Eq. ([53]) can be calculated indepen- 
dently with the coupling parameter approach described 
in Sec. [IVJ This provides a precise control of any ap- 
proximation involving the cutoff /3q. We did compute 
Fq (to, L = a (L z — |) , = aL x , S+ = as+J within the 



12 



aforementioned coupling parameter approach and we 
have taken into account this term in Eq. (|53p whenever 
it is relevant within the statistical precision. The nu- 
merical integrations in Eq. (1541) have been carried out 
according to Simpson's rule. Certain technical details 
are reported in App. O Finally, the determination of the 
critical Casimir force on the basis of Eq. ((53)) requires 
the knowledge of the reduced bulk free energy density 
-Fbuik(/3) which is independent of the BC. We have de- 
termined it via MC simulations of lattices size L\ with 
L z = 24 — 256 and periodic BC. In App. [B] we report 
certain details of this computation, which is important 
for a successful determination of Fq ■ 

Along these lines we have computed the critical 
Casimir force for lattice thickness L z — 8, 12, 16, and 
24 with the BC shown in Figs. [1] and [3] as well as for 
k = 0, 1/2, 1, 2, and 3. As in Sec. |IV] we have consid- 
ered three aspect ratios for each value of L z and n: we 
have taken p = 1/8, 1/12, and 1/16 for k < 2, as well as 
p = 1/12, 1/18, and 1/24 for n = 3. We have checked 
that for these small values the data are independent of p 
within the statistical accuracy. Therefore we expect that 
our results capture the limit p — > 

In the present case, for t ^ it is not easy to subtract 
the scaling corrections because according to Eq. (fT9|) a 
part of the scaling corrections oc 1/ L z stem from the de- 
pendence on L z of the second scaling argument of 9. This 
holds even if the scaling ansatz of Eq. (|19p does not com- 
pletely capture the 1/L Z scaling corrections: in fact, the 
non- universal length c ', defined in Eq. (|33[) and extracted 
from the fits reported in Tables U and [TT1 shows a small 
but significant dependence on k, which would be absent 
if scaling corrections were independent of k. In Ref. [2(| a 
similar problem was encountered in the MC investigation 
of the critical Casimir force in the presence of an isolated 
chemical step. There the dependence of the force on the 
aspect ratio contributes to the scaling corrections. Since 
this dependence on p was found to be linear, in that case 
it was possible to eliminate the scaling corrections via a 
first-order Taylor expansion of the critical Casimir force 
in p. As Figs.[5]and|6]show, in the present case the critical 
Casimir force does not follow such a simple dependence 
on k. Furthermore, the possible values of n which can be 
sampled by the MC simulations are constrained by the 
fact that the stripe width s + has to be an integer number. 
Due to these technical difficulties, here we implement an 
approximate scheme for the removal of the scaling correc- 
tions. For every value of k we extract the non- universal 
length c ' from the fits of Tables U and HU Then we employ 
the substitution L z — > L z + c ' . Since such a substitu- 
tion cannot completely eliminate the scaling corrections 
oc L z l , the resulting scaling function 6{t,k) exhibits a 
residual scaling correction oc ip{r, k)/L z , where ip(r, k) is 
a scaling function. By construction, we have ip(0, k) = 0. 
Thus, since ip is a continuous function, there is an interval 
around r — in which the residual scaling corrections are 
negligible with respect to the numerical precision. Fur- 
thermore, for k — t and k — > oo this method becomes 




T 

FIG. 7. The universal scaling function #( + d )(t) of the critical 
Casimir force for the BC (+, o) shown in Fig.0 corresponding 
to the limit k = s+j L z — > of the BC shown in Fig.Q] Scaling 
corrections have been subtracted by using c ' = 0.36(4) (see 
the main text). We also compare our results with those of 
Ref. [H for L = 16 and of Ref. [H for L = 20. The omitted 
statistical error bars are, apart from r < —10, comparable 
with the symbol size. 



exact and thus we have ip(r, k 0) = ip(r, k — > oo) = 0. 
Therefore, the interval of validity around r = is ex- 
pected to increase as n is lowered towards or is increased 
toward oo. 

In Fig. [7] we show our results for the BC shown in 
Fig. |31 corresponding to the limit n = s + /L z — > of the 
BC shown in Fig. [1] In order to normalize the scaling 
variable r, one needs the value of the non-universal am- 
plitude £g] of the correlation length £. From Ref. [2j| we 
infer £^ = 0.4145(4) in units of the lattice constant. As 
for the critical exponent v, we use the recent MC result 
v = 0.63002(10) of Ref. In FigJT] we also com- 

pare our results with those of Refs. [31( and [32]]. We 
observe a perfect agreement with the results of Ref. [3l[ , 
which in fact have been obtained by simulating precisely 
the same improved Blume-Capel model. The compari- 
son with the results of Ref. 32] is less satisfactory and 
reveals a difference between the curves around the po- 
sition of their maximum in the low-temperature phase, 
i.e., r < 0. This difference may be due to the fact that 
the Ising model simulated in Ref. [32] suffers from larger 
scaling corrections than the improved model used here, 
which makes the extrapolation of the FSS limit more dif- 
ficult. For the BC illustrated in Fig. [Q in Figs. [81 151 [101 
[Til and we show our results for the scaling function 
6+(t,k), for k = 1/4, 1/2, 1, 2, and 3, respectively. 

Inspection of Figs. I71TT21 reveals a satisfactory scaling 
collapse for the lattice sizes considered here. This sup- 
ports the validity of the procedure described above to 
eliminate the scaling corrections. In Fig. \T3\ we show a 
comparison of the critical Casimir force for k = 0, 1/4, 
1/2, 1, 2, and 3, as obtained for L z — 24. We also com- 
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FIG. 8. The universal scaling function 0+(t, k) of the critical 
Casimir force for the BC shown in Fig. Q] with k = S+ /L 
I 1 and c ' = —0.48(2). The omitted statistical error bars are 
comparable with the symbol size. 



FIG. 11. Same as Fig.[S]for k = 2 and c 



-0.14(1) 





FIG. 12. Same as Fig.fjfor k = 3 and c ' = -0.10(1). 



FIG. 9. Same as Fig.rjfor n = 1/2 and c ' = -0.44(1). 




FIG. 10. Same as Fig.fjfor k = 1 and c ' = -0.26(1). 



pare the present results with the universal scaling func- 
tion which describes the critical Casimir force for an iso- 
lated chemical step in the limit of vanishing aspect ratio, 
as determined in Ref. (20j . This system corresponds to 
the limit k — > oo and results in the mean value of the crit- 
ical Casimir force for laterally homogeneous (+, +) and 
(+,— ) BC. In the whole range < n < oo the critical 
Casimir force is always repulsive. This is expected be- 
cause the stripe width for (+) and for (— ) BC are equal 
and the repulsive critical Casimir force for (+, — ) BC 
is stronger than the attractive one for (+, +) BC - 
In Fig. [T3] we also show a comparison with the mean 
value of the critical Casimir force for the homogeneous 
(+, +) and [+, — ) BC, as obtained by MC simulations in 
Refs. [HHl. 

In Fig. [14] we show our results for the (o, o) BC shown 
in Fig. [H corresponding to the limit k = s+/L z — > of 
the BC (o) vs stripes shown in Fig. [5] We also compare 
our results with those of Ref. [25[ for the approximants (i) 
and (ii) presented therein. The approximant (i) agrees 
with our results for r > — 6, whereas the approximant 



14 




10 15 



-0.5 



K 

o 

® -1.5 



-2.5 




r 



V 



* tit'* 



""-20 



-5 
X 



(o) vs stripes 



• L z = 


8,p = 


1/8 


■ L z = 


8,p = 


1 / 12 


* L z = 


8,p = 


1 / 16 


* L z = 


12, p 


= 1/8 


- L z = 


12, p 


= 1/12 


L 7 = 


12, p 


= 1/16 


► L,= 


16, p 


= 1/8 


L z = 


16, p 


= 1/12 


" L z = 


16, p 


= 1/16 



10 



FIG. 13. Comparison of the universal scaling function 
8+(t,k) for k = 0, 1/4, 1/2, 1, 2, and 3 as determined with 
L — 24. We compare the data also with the scaling func- 
tion 8+ (t, k — > oo ) in the limit of vanishing aspect ratio p, 
as obtained in Ref. 20] with L — 16. The limit k — > oo 
corresponds to the critical Casimir force between a homoge- 
neous (+) surface and a surface with an isolated chemical 
step which, for p — > 0, results in the mean ualue of the criti- 
cal Casimir force for laterally homogeneous (+, +) and (+, — ) 
BC [20(|. We compare the results also with those latter mean 
values, which are either extracted from the so-called approxi- 
mant IV of Ref. [25[ (mv (IV)) or which stem from the results 
of Ref. @ (mv). 
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FIG. 14. Universal scaling function 0( O]O )(t) of the critical 
Casimir force for the BC (o, o) shown in Fig. [4] corresponding 
to the limit k = s+/L z — > of the BC shown in Fig. [2] 
Scaling corrections have been subtracted by using c ' = 0.8(2) 
(see the main text). We compare our results with those of 
Ref. [2^] obtained from the approximants (i) and (ii) presented 
therein and for the film thickness L — 20. The inset provides 
a magnification of the resulting curves close to the minimum 
of the force, for the largest available film thickness L — 16 
and for the three aspect ratios p = L z /L x considered here. 



FIG. 15. Universal scaling function o (t,k) of the critical 
Casimir force for the BC shown in Fig. [5] with k = s+/L z = 
1/4 and c ' = 0.02(9). The data points for L = 8 and p = 
1/8, 1/12 are hardly visible because they overlap with the 
other data sets. The inset provides a magnification of the 
resulting curves close to the minimum of the force, for the 
largest available film thickness L — 16 and for the three aspect 
ratios p = L z /L x considered here. 





-0.5 
-1 

ii -1.5 

o 

-2 
-2.5 
-3 



(o) vs stripes 



K=l/2 



-2.5 



'-20 



1 » ' 

-15 



• L z = 


8,p = 


1/8 


■ L z = 


8,p = 


1 / 12 


L z = 


8,p = 


1 / 16 


» L z = 


12, p 


= 1/8 


- h z = 


12, p 


= 1/12 


L z = 


12, p 


= 1/16 


' L z = 


16, p 


= 1/8 


L z = 


16, p 


= 1/12 


" L z = 


16, p 


= 1/16 



-10 



-5 
X 



10 



FIG. 16. Same as Fig. Hoi for k = 1/2 and c ' = 0.05(8). 



(ii) displays a systematic deviation from our results. For 
t < — 6 both approximants show a disagreement with 
our results. While the approximant (ii) displays a small 
but visible deviation from our results, the approximant 
(i) exhibits a larger, systematic deviation from our re- 
sults. Such deviations may be due to the difficulty in 
extrapolating the FSS limit of the Ising model used in 
Ref. [25|, which exhibits larger scaling corrections than 
the improved model of Eq. ([To]) . For the BC illustrated 

in Fig. [I in Figs. EI OS EH [H and ^ we show our 
results for the scaling function o (r, k), for k = 1/4, 1/2, 
3/4, 1, and 3, respectively. 

The numerical determination of the critical Casimir 
forces in the presence of a Dirichlet BC at one of the two 
confining surfaces has turned out to be much more in- 
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FIG. 17. Same as Fig. [15] for k = 3/4 and c ' = 0.37(7). 
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FIG. 18. Same as Fig.[To]for « = 1 and c ' = 1.3(2). 
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FIG. 19. Same as Fig. [15] for k = 3 and c ' = 0.36(9). We 
also compare our results with the interface estimate given by 
the r.h.s. of Eq. (|15p . The scaling function changes sign at 
T = TO ~ -2.7. 



volved than the computation for the BC of Figs. Q] and [3] 
First, at variance with the previous cases, we observed 
the onset of a dependence of the critical Casimir force on 
the aspect ratio p = L z j L x . As illustrated in the insets 
of Figs. [F4"lll9[ such a dependence on p appears in a nar- 
row interval of r in the low-temperature phase. Although 
small, the differences between the calculated scaling func- 
tions o (t, k) for the three aspect ratios considered here is 
visible and larger than the statistical error bars [6^ . The 
observed dependence on p implies the onset of a lateral 
correlation length, associated with an ordering process 
in the low temperature phase. In order to understand 
this point, it is useful to consider the limit /? — > oo, i.e., 
the ground state of the model with the BC illustrated in 
Figs. [2] and [4] For the BC shown in Fig. [4] it is easy 
to see that the ground state is a spatially homogeneous 
state in which all spins take the same value. For the BC 
shown in [21 besides the homogeneous state, one can con- 
sider also a "striped" state, in which each spin in the film 
takes the value corresponding to the underlying stripe, so 
that the configuration of the system consists of columns 
of cross sectional area s+ x L x and height L z . In view of 
the periodic BC in the two lateral directions, the area A 
of the interface between + and — spins is given by 



A = 



2 



L L 2 
A = — L Z L X = — j 



homogeneous state, 
striped state. 



(55) 



Thus, at low temperature, the system orders in a homoge- 
neous state for n < 2 and in a striped state for n > 2. As 
a function of the parameter k, the ground state undergoes 
a first-order transition at K = 2. Moreover, for k = 2, 
besides the homogeneous and the striped ground states, 
there are other states which have the same (minimal) en- 
ergy: such states can be obtained by flipping the value 
of the spins in a single column in the striped state. We 
note that the number of these additional ground states 
diverges in the thermodynamic limit. The emergence of 
these ground states at k = 2 gives rise to a sort of glassy 
behavior at low temperatures, which results in a con- 
siderable technical difficulty in simulating these systems. 
We leave this issue for future research. 

This lateral ordering process at low temperatures cor- 
responds to a phase transition which occurs in the film ge- 
ometry characterized by the BC described by Figs.[5]and 
2] This causes the dependence on the aspect ratio exhib- 
ited in Figs. [2][T!J] We note that, for the BC correspond- 
ing to Figs. [1] and [3[ the striped state is never a ground 
state. Moreover, without an external bulk field the pres- 
ence of a surface field at the upper surface rounds the 
transition between the paramagnetic high-temperature 
phase and the homogeneous ground state to a simple 
crossover. This is in agreement with the independence 
of p observed in Figs. [7][T2l The appearance of a lateral 
correlation length breaks the scaling behavior discussed 
in Sec. [II] On the other hand, inspection of Figs. [T4l 
[12] reveals that the data for the two smallest aspect ra- 



tios agree within the statistical error. Therefore, since 
one expects a smooth dependence of the scaling function 
9 (t,k) on p, in particular in the limit of p — > 0, we 
can regard our results for the smallest aspect ratio as a 
reliable extrapolation of the limit p — > 0. 

Another difficulty in the numerical determination of 
the critical Casimir force for the BC shown in Figs. [5] 
and S] lies in the fact that the scaling function 9 (t, k) ex- 
hibits a minimum in the low-temperature phase which is 
shifted towards more negative values of r upon increasing 
k. Thus, in order to study this important feature of the 
scaling function, one has to generate MC data for tem- 
peratures lower than the ones needed for the BC shown 
in Figs. [T] and O Upon lowering the temperature the sim- 
ulations become increasingly difficult because of the ap- 
pearance of many metastable states associated with the 
aforementioned ground-state phase transition at K = 2. 

Finally, in order to eliminate the leading scaling cor- 
rections, we have implemented the procedure outlined 
above. We note that for the BC shown in Fig. [2] such a 
method appears to be less reliable. While for k < 1/2 
and k — 3 the overall scaling collapse is good, for n = 3/4 
and for sufficiently negative values of r, there is a small 
but systematic deviation between the data for lattice size 
L = 12 and L = 16. The scaling collapse is even worse 
for k = 1; in this case a further complication seems to 
be that, apparently, in this case scaling corrections are 
stronger (see Table ITV| . 

According to the discussion in Sec. IIII1 for the BC 
shown in Fig. [5] in the limit n — > oo one expects to re- 
cover the BC shown in Fig. [3] Since for k — the force is 
always attractive (see Fig. I14|) and for k —> oo the force is 
repulsive (see Fig. [7]), at a certain intermediate value of n 
the force has to change sign. According to Fig. [BJ at crit- 
icality this occurs at k — kq ss 1.2. Besides a change of 
sign of the force as a function of k there is also a change 
of sign as a function of r. This is nicely illustrated in 
Fig. Q32 where for k = 3 the force is found to be repul- 
sive (resp. attractive) for r > tq (rcsp. r < To), with 
To ~ —2.7. This implies that in the scaling regime and 
for a given temperature T < T c , i.e., t = (T-T c )/T c < 
the force is repulsive (resp. attractive) for L < L (t) 
(resp. L > L (t)), with L (t) = ^oi(T /t)". Therefore 
L = Lo(t) is a mechanically stable point of equilibrium 
for the critical Casimir force which can be sensitively 
tuned by varying the reduced temperature. This can be 
exploited for levitation purposes [41[ . In Fig. Q1J] we also 
compare our result with the interface estimate, i.e., the 
r.h.s. of Eq. (|15l) . which is expected to hold for k > 2 
and r <C — 1. To this end, we employ the estimate of the 
universal amplitude ratio R a — 0.377(11) (65|. The in- 
terface estimate is in nice agreement with our MC results 
for r < -3.5. 

In Fig. [30] we show a comparison of the scaling function 
o (t, k) of the critical Casimir force for the BC shown in 
Fig. [2] for k = 0, 1/4, 1/2, 3/4, 1, and 3 as determined 
with L = 16 and with the smallest aspect ratio p avail- 
able. We also compare these results with the Casimir 
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FIG. 20. Comparison of the universal scaling function 
9 {t,k) for k = 0, 1/4, 1/2, 3/4, 1, and 3 for the BC (o) 
vs stripes shown in Fig. [2] as determined with L = 16 and 
the smallest aspect ratio p available. We compare these data 
also with the scaling function 9 (t,k — > oo) = #( +j0 )(t), as 
obtained in Fig. [7] with L = 24. For further discussions see 
the main text. 

scaling functions for the BC (+, o) shown in Fig.[3l which 
corresponds to the limit k —> oo. Figure I2TJ1 suggests that 
the approach of the limit k — > oo is somehow singular. 
Apparently, for every finite value of k, the force becomes 
attractive for sufficiently negative values of r and exhibits 
a minimum which deepens and shifts to more negative 
values of t as k is increased. Simultaneously, the zero of 
o (t, k) shifts towards lower values of r. 

VI. MEAN-FIELD THEORY 

Within the field-theoretic approach, bulk and surface 
critical phenomena of the Ising universality class are de- 
scribed by the standard Landau- Ginzburg- Wilson fixed- 
point Hamiltonian given by d, |66[ 

/^d^Vj^-M^}, (56) 

where </>(r) is the spatially varying order parameter 
describing the critical medium, which completely fills 
the volume V bounded by the boundaries dV in d- 
dimensional space. In Eq. (|56"|) f oc t and u > is a 
coupling constant providing stability for t < 0; c(r) is 
the surface enhancement, which, within mean-field the- 
ory, can be interpreted as an inverse extrapolation length 
of the order parameter field, and h\(r) is an (external) 
surface field acting on the order parameter at the bound- 
aries. Here, we consider surface fields and enhancements 
which can be different for the two confining surfaces and 
which may also vary along one lateral direction of a single 
surface. In the strong adsorption limit, i.e., (±) BC, cor- 
responding to the so-called normal surface UC, the sur- 
face behavior is described by the renormalization-group 
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fixed-point values hi —> ±00, and the order parameter 
diverges close to the surface: <f)\av — > ±00. The or- 
dinary surface UC corresponds to the fixed point val- 
ues {c = 00, hi — 0} and a vanishing order parameter 
4>\dv — 0, i.e., Dirichlet (o) BC. The film geometry con- 
sidered here is bounded by surfaces at z = and at z = L 
with either homogeneous (+) or (o) BC or periodically 
alternating (+)/(-) BC of width S+ = P/2 along the 
lateral z-direction (see Figs.[TH2J). 

The Hamiltonian given in Eq. (|56l) is minimized by 
the mean-field order parameter profile m = u 1 / 2 (</>): 
SH[(f)]/&(f)\ r j )=u -i/2 m = 0. Renormalization group argu- 
ments tell that mean-field theory (MFT) provides the 
correct universal properties of critical phenomena for 
spatial dimensions above the upper critical dimension 
d > due = 4 (up to logarithmic corrections in d — d uc ). 
Mean-field theory provides the lowest order contribution 
to universal properties within an expansion in terms of 
4 — d = e. Thus, universal properties in d = 4 can 
be determined from MFT, up to two independent non- 
universal amplitudes appearing in the description of bulk 
critical phenomena (two-scale universality [J, |5|): the 
amplitude B of the bulk order parameter ((f)) = ±_B|i|' 3 
for t < 0, where /3(d = 4) = 1/2, and the amplitude £q~ 
of the correlation length [see Eq. ©, where v{d = 4) = 
1/2]. Since here we are dealing only with vanishing or 
diverging values of hi and c, within MFT all quantities 
appearing in Eq. (l56l) can be expressed in terms of these 
amplitudes: f = t($^)~ 2 and u = 6(B£ i q)^ 2 . Using the 
stress tensor method [2l| the mean-field universal scaling 
functions of the critical Casimir forces at the upper crit- 
ical dimension d uc = 4 can be inferred directly from the 
MFT order parameter profiles up to an overall prefactor 
oc it -1 . 

For the laterally homogeneous (+,+), (+,—), (+,0), 
or (o, o) BC the MFT order parameter profiles across 
the film [U, [(J]} and the corresponding universal scaling 
functions of the critical Casimir force are known analyt- 
ically [2]], [68[ . Accordingly, the critical Casimir ampli- 
tude 9 (+)+) = 8if 4 (l/v / 2)(S^+) 2 ~ -47.2682(B£+) 2 , 
where K(k) is the complete elliptic integral of the first 
kind 12.11. Note that, within MFT, the scaling func- 



tions (+ ,_)(t) = -40 (+ , +) (-t/2) [69|, and (+ , o) (r) = 
0r + _)(4r)/16 [21j are directly related to each other, so 



(+,0) 



-0 



(+,+) 



/4- 



7(_)_ {^t ) / iu are uuecuy reiaiei 

that at T c 0( +) _) = -40 (++ ) and 6 
In contrast to the case d = 3, the MFT scaling function 
for (0,0) BC vanishes for r > (i.e., 9( 0j0 )(d = 4) = 0) 
and exhibits a cusp-like singularity at its minimum at 
t = — 7r 2 below which 0( o . o )(t < —n 2 ) = #(+.+) (t) and 
above which an analytic expression for #( , ) has been 
derived in Ref. (r38|. 

In order to obtain the spatially inhomogeneous MFT 
order parameter profile for the film geometry involving 
chemically striped surfaces, we have minimized H[4>] nu- 
merically using a quadratic finite clement method. Here, 
we extend previous investigations [38| to negative values 
t < and to a broader range of geometrical parame- 
ters. The corresponding scaling functions for the critical 
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FIG. 21. Reduced critical Casimir force amplitude 0+(k) 
[Eq. JTTJ] in units of |6( +j+) | for the BC shown in Fig. Q] 
as obtained within mean-field theory. For k, —¥ the Casimir 
amplitude approaches the value for (+, o) BC shown in Fig. [3] 
i.e., O( + , o )/|0( +i+ )| = j, indicated by the lower red dotted 
line. For large stripes, 6+(k — > oo)/|0( + + )| approaches the 
average value of the reduced Casimir amplitudes for (+, +) 
and (+,-) BC, i.e., (0 (+i+) + e (+i _ ) )/|20 (+ , +) | = § shown 
as upper blue dotted line. For « 3> 1 the behavior of the 
Casimir amplitude 0+(/t)/[0( +i+ )| approaches the function 
I — I a + K _1 (see the black dashed line and the main text). 
From a least square fit we have obtained a, — 0.420(4). Com- 



pare Fig. where 0( + )/ 
oo)/|0 (+ , +) |~ 2.91(5). ' 



|0(+,+)| = 0.60(1) and Q+{k 



2il. 



Casimir force are obtained via the stress tensor 

The boundary condition for the diverging order param- 
eter profile at those parts of the surface where there are 
(+) or (— ) BC can be implemented numerically only ap- 
proximately via a short-distance expansion of the corre- 
sponding profile for the semi-infinite systems 0,0 ■ Thus, 
the MFT data presented below are subject to a numeri- 
cal error which contains also the uncertainties due to the 
fineness of the numerical mesh. We estimate the numeri- 
cal error for the data presented below to be less than 1% 



or ±0.004 x 16 



(+. 



if the latter is bigger. 



A. Critical Casimir amplitude at T c 

In Fig. [5T] the amplitude of the critical Casimir force 
9+(k) = + (0,k) (see Eqs. © and (fTTjll for a striped 
surface opposite to a homogeneous surface with (+) BC 
is shown as obtained numerically within MFT in units 
of |0( +i+ )|. We have been able to calculate the val- 
ues of 0-|_ (k) numerically within the range k = 0.1 to 
k = 80. As discussed above, for k —> the Casimir 
amplitude approaches the value for (+,0) BC shown in 
Fig. El i.e., 0(+, o ); so that for relatively narrow stripes 
the chemically striped wall effectively mimics a wall with 
(o) BC. On the other hand, for n — > 00 the Casimir 
amplitude approaches the average value of the Casimir 
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FIG. 22. Reduced Casimir amplitude 9 (re) [Eq. ([TT])] in 
units of |0( +i+ ) | for the BC shown in Fig.[2]as obtained within 
MFT. For re —¥ the Casimir amplitude approaches mono- 
tonically from positive values the limiting value 0( o , o ) = 
shown by the lower green dotted line. According to Eq. (|14|1 , 
for re ^> 1 the reduced Casimir amplitude o (re)/|0(+,+)| ap- 
proaches |(1 — q k _1 ) shown as black dashed line. From a 
least square fit we have obtained, within MFT, a = 0.857(9) 
[Eq. ((58)) ]. For re — > oo, o (re)/|0( +i _|_)| approaches the 



Casimir amplitude for (+,o) BC, i.e., 



(+,o)/ 1©(+,+) 



1/4 



shown as the upper red dotted line. Compare Fig. [5] where 
o (re = 0)/|© (+j+) | = 0.037(6) and o (re -> oo)/|0 (+ , +) | = 
0.60(1). 



amplitudes for (- 



and (+, — ) BC, i.e., 0+( K ~^ °°) = 



(e (+ , +) + e (+ ,_))/2 = -fe 



(+,+)) whereas _)_(«) mono- 
tonically interpolates between these two limits. 

For /?> 1, according to Eq. (fl"3|) . we find for the crit- 
ical Casimir amplitude 



0+(re> 1) 
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2 4 + 



(57) 



where the proportionality constant a + has been deter- 
mined within MFT as a + = 0.420(4) by using a least- 
square fit. 

Figure [521 shows the reduced critical Casimir force 
amplitude o (re) in units of |0( + . + )| for the case of a 
striped surface opposite to a surface with a homogeneous 
(o) BC (see Figs. and gj). Similar to Fig. ED ©o(«0 
monotonically interpolates between the limiting values 



for k 



and re 



oo, i.e. 



e 



and 



0( +iO )/|0( +i+ )| = 1/4, respectively. For small stripes 
the amplitude o (k — > 0) approaches its limit already 
for larger values of re than in the case of a homogeneous 
(+) BC shown in Fig.[5T] This indicates that the strength 
of the tendency of a chemically striped surface to effec- 
tively mimic an (o) BC in the limit re — > also depends 
on the type of homogeneous BC at the opposing surface 



FIG. 23. Reduced universal scaling function 

0+(t, /t)/|0(_|_ i+ )| [Eq. @] for a striped surface opposite to 
a surface with homogeneous (+) BC (Fig. [TJ, as deter- 
mined numerically within MFT for various values of re. For 
re — v and re — > oo, the reduced scaling functions ap- 
proach their limiting behaviors #(+,<,) ( r )/|0(+.+) i [Eq. l|12fl] 
and (*(+,+) (t) + #(+,-) (r))/|2e (+ ,+) | [Eq. respectively. 
Compare Fig. [13] by taking into account that there, i.e., in 
d = 3, |G (+>+) | = 0.820(15). 



of the film. According to Eq. ([14"|). for re 3> 1 the depen- 
dence of the Casimir amplitude o (re) on re approaches 
the following form: 

e a ( K » i) ~ (o , o) + (0 (+jO) - © (0)0) ) (i - ^) 







(+,+) 



(1 - a k- 1 ) , (58) 



where we have determined a a = 0.857(9) via a least- 
square fit. 

Whereas the behavior of the Casimir amplitude 0+ (re) 
for the case of a homogeneous (+) BC as calculated 
within MFT fFig. l2"Tj) is similar to the one obtained from 
MC simulations (Fig. [5]), the form of O («) for the case of 
a homogeneous (o) BC as obtained within MFT (Fig. |2"2")) 
is qualitatively different from the one obtained from MC 
simulations (Fig. [6]). This will be addressed in more de- 
tail in Sec. IVIII below. 



B. Scaling function of the critical Casimir force 

The reduced scaling function 6+{t, /t)/|0(+ ) +)| 
[Eq. ©] of the critical Casimir force between a chem- 
ically striped surface and a homogeneous surface with 
(+) BC (Fig. HJ is shown in Fig. OH for d = 4 (MFT) 
and for various values of re. For re — > 0, 9+(t, k)/|0(+.+)| 
approaches the scaling function 6*(+, o )(' r )/|0(+,+)l> i- e -> 
the striped surface effectively mimics a surface with 
homogeneous (o) BC. On the other hand, for re — > oo, 
the universal scaling function of the critical Casimir 
force approaches the average of the scaling functions for 
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FIG. 24. Reduced universal scaling function 

9 (t, k)/|0/_|_ + n| of the critical Casimir force for a striped sur- 
face opposite to a surface with homogeneous (o) BC (Fig. [2J, 
as determined numerically within MFT for various values of 
re. We compare the data also with the reduced scaling func- 
tions 0( O]O )(t)/|0( +>+ )| and 0( +iO )(t)/|0( +i+ )|, which corre- 
spond to the limits re -> [Eq. (fT5|) ] and re — s> 00 [Eq. ((14)) ] . 
respectively. For re < 2 the numerically obtained MFT scaling 
functions suggest the occurrence of a cusp- like singularity or a 
finite jump of 9 (r, re) at its minimum position r m i n . (Due to 
the numerical difficulties in determining the thermodynami- 
cally stable configuration, both the positions and the depths 
of the minima of o (t,k < 2)/[0( +i+ )| are affected by an es- 
timated numerical error of around 10%, which is one order 
of magnitude larger than for the remaining data.) For re > 2 
the scaling functions diverge for r —5- —00 [Eq. (| 15f) ] . Com- 
pare Fig. [20] by taking into account that there, i.e., in d — 3, 
|0 (+>+) |= 0.820(15). 



(+,+) and (+,-) BC, i.e., 0+(t,k -> oo)/|0 (+)+) | = 
(0 (+ , +) (T)+0 (+ ,_)(T))/|2e (t)+) | [Eq. (EM. For inter- 
mediate values of n, the scaling functions smoothly and 
monotonically interpolate between these limiting cases. 

As discussed in Sec. [VI the behavior of the universal 
scaling scaling function 8 (r, re) for a striped surface op- 
posite to a surface with homogeneous (o) BC (Fig. ^ is 
more complex than the one in the previous case. Whereas 
for t > the scaling function 8 q (t, k) smoothly interpo- 
lates between its limiting behaviors 0( o ,o){ T ) f° r K = 
and #(+.o)( T ) for K °Oi for negative values of r its 
dependence on re is non-monotonic and involves a phase 
transition associated with the one at re = 2 between the 
ground states of the system (see Eq. ([55)) ). For re < 2 the 
ground states are spatially homogeneous, which results in 
a vanishing value 8 (t — > —00, re < 2) — > 0. The numer- 
ically obtained MFT data shown in Fig. [24] suggest that 
the minima of the scaling functions for re < 2 correspond 
to a cusp-like singularity or even a finite jump. (Re- 
call that 8 is the scaling function of the critical Casimir 
force, which is the derivative of the Casimir interaction.) 
However, due to the presence of metastable striped and 
homogeneous states the numerics even within MFT is so 
involved that the present data suffer from an error of the 
position of the minimum of around 10%. Moreover, due 
to using the short-distance expansion in the numerical 



implementation of (±) BC, it is technically difficult to 
distinguish these metastable states for re ~ 2. For n > 2 
a striped ground state is stable, which involves a diver- 
gence of the scaling function for t — > —00 so that for 
t < the transition to its limiting behavior 6*( + q )(t) > 
for k — y 00 is somewhat singular. Since at T — T c , the 
critical Casimir amplitude o (k) is non- negative for all 
values of k (see Fig. [5U for re < 0.5, O is vanishingly 
small), within MFT the scaling function 8 (t, k) changes 
sign for all values of re at a certain value t*(k) < 0. 

In the following we consider the contribution of the 
interface tension to the critical Casimir force for t < 
(see Eq. (|15[) ). Near T c the interface tension varies as 
c = <7o|£| M where fj, = (d — so that /i = 3/2 
within MFT [65| : o~q is the corresponding non-universal 
amplitude which forms the universal amplitude ratio 
TGT^oKo")^ = R*- Within MFT a/(k B T c ) = 
^ti- 1 ^)-^ 1 )!*]" so that R a = | N /2(B^+) 2 and 
i2 ff /|0(+ ) +)| — 0.020. For the homogeneous configuration 
with the interfaces parallel to the film (i.e., for re < 2), 
the interface energy does not contribute explicitly to the 
resulting force because the area of these interfaces is not 
changed upon varying of the film thickness. (Note, how- 
ever, that the order parameter profile across these inter- 
faces does depend on L.) For the striped configurations, 
i.e., for k > 2, in which the interfaces are oriented perpen- 
dicular to the film, the interface tension dominates the 
resulting force for large negative r (i.e., L large), because 
approximately the interface along the z-direction has an 
area L^~ 2 L which is proportional to the film thickness 

L. Thus, the free energy T\ of such a single interface is 
given by 

= L^Lcr, (59) 

where Lu is the extension of the system along the invari- 
ant direction(s). For a single such interface this gives rise 
to a force along the normal direction: 

Fi s = ~± = -L\- 2 o~. (60) 

For the striped state there are 2 x L\\/P = Lh/S+ such 
interfaces so that the total force per area Ly 1 of the film 
and per ksT c is 

F T,tot = 1 O" = 1 1 rd-X °0 U\(d-l)v 

k B T c VV x S+k B T c L d n k B T c ] 1 

(61) 

so that its contribution 8 q j-{t, re) to the universal scaling 
function of the critical Casimir force reads (see Eq. ([9])) 

B , t {t,k) = -^\tT, (62) 
re 

which is attractive and becomes as strong as 0( + . + j for 
|r| M /re > 50 within MFT. Accordingly, for the limit 
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FIG. 26. Comparison of the normalized critical Casimir 
amplitude &+(k) = [0+(k) - <3>+{k -> 0)]/[6+(k ->• oo) - 
0+(k — >• 0)] [Eq. (16311 ] for a homogeneous (+) wall opposite 
to a striped wall (Fig. [TJ as obtained from MC data (symbols; 
same as Fig. [5]) and within MFT (solid line; see Fig. 121]). 



VII. COMPARISON BETWEEN MEAN-FIELD 
THEORY AND MONTE CARLO DATA 



FIG. 25. Reduced universal scaling function 

9 (t, k;)/|0( + . + ) | of the critical Casimir force for a striped sur- 
face opposite to a surface with homogeneous (o) BC (Fig. [2]), 
as determined numerically within MFT (solid lines, same as 
Fig. |2~4")| . For r <C — 1 and k > 2 they agree well with 
the asymptotic expression given by the r.h.s. of Eq. (|15[) 
shown as dashed lines (a). For k > 2 and large negative 
values of r, i.e., r <C —10, the attractive interface contribu- 
tion — i? CT K _1 |r| M /|0( + , +) | [Eq. ([62l ] dominates the the scal- 
ing function o (t, k)/|6( +j+ )| (b). 



t <C — 1 and k > 2 the scaling function of the crit- 
ical Casimir force approaches the expression given in 
Eq. (fl~5)) . which corresponds to the sum of the homo- 
geneous contribution and the contribution due to the in- 
terfaces oriented perpendicular to the film surfaces. 

Figure [231 compares 8 (t, k) for a striped surface oppo- 
site to a surface with homogeneous (o) BC as determined 
numerically within MFT with the estimate of the corre- 
sponding interface contribution as given in Eq. (|T5|) . The 
dashed lines shown in Fig. [55] correspond to Eq. (fTol) . 
They are approached by the actual scaling functions 
shown as solid lines in Fig. [551 As expected, Eq. (fT5|) 
does neither describe the behavior for n < 2 nor the one 
for small absolute values of r. However, for r -C — 1 and 
K > 2, the scaling functions agree rather well with their 
asymptotic behavior given in Eq. (]15[) . 



A. Critical Casimir amplitude at T c 

Different than the MC data for d — 3, the universal 
scaling functions of the critical Casimir force obtained 
within mean-field theory can be determined only up to 
an unknown constant amplitude. In order to facilitate 
nonetheless a valuable comparison between them, which 
illustrates the dependence of the scaling functions on the 
spatial dimension d, it is useful to normalize them by 
an overall amplitude so that the unknown constant am- 
plitude for the MFT results drops out. In the previous 
section we normalized the various scaling functions by 
one and the same universal critical Casimir amplitude 
|G( +j _|_)|. Here, we propose an alternative normalization, 
which makes use only of that scaling function under con- 
sideration and also normalizes the ratios between the cor- 
responding critical Casimir amplitudes, which depend on 
d: 

a, . e(K)-e(«-»o) r o, k^o, 

9(k ^ oo) - 9(k ^ 0) \l, k^oo. 

(63) 

As discussed in the previous sections, the critical 
Casimir amplitude between a chemically striped wall 
and a homogeneous wall with (+) BC interpolates be- 
tween 0+(k — > 0) = ©(+, ) an d - > o°) = 
(0( + + ) + 0( +i _))/2. Figure [5FH shows the corresponding 

normalized critical Casimir amplitude 0+(k) [Eq. (|63|) ] 
as obtained from MC data (symbols) as well as obtained 
within MFT (full line). As can be inferred from Fig. [231 
the behavior of the normalized Casimir amplitude + (k) 
as a function of k as obtained from MFT (d = 4) is rather 
similar to the one in d — 3. Thus, for this geometry the 
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FIG. 27. Normalized Casimir amplitude o (/-c) = [0 o (k) — 
9 (k -> 0)]/[6 o (k -> oo) - 9„(k ->■ 0)] [Eq. JB3J] for a ho- 
mogeneous (o) wall opposite to a striped wall (Fig. [2]) as ob- 
tained from MC data (symbols; same as Fig. [S]) and within 
MFT (solid line; see Fig. I22[) . In contrast to the behavior 
shown in Fig. 1261 the MFT results differ qualitatively from 
the behavior in d = 3. In both cases MFT overestimates the 
strength of the force (here for k > 0.75). Q (k — > oo) attains 
its limiting value 1 slower than 0+(k — > oo). 



effects of the chemical patterning are captured even semi- 
quantitatively by MFT. 

In contrast, for the case of a homogeneous (o) surface 
opposite to a striped one (Fig. [2]), we find qualitative dif- 
ferences. In Fig. 1271 the normalized critical Casimir am- 
plitude o (k) [Eq. (1S3"]) ]. as obtained both in d = 3 and 
within MFT, is shown, using the corresponding limits 
& {k —> 0) = O( ,o) and 6„(k —> oo) = 0(+, o ). Whereas 
the critical Casimir amplitude as obtained from MC sim- 
ulations shows a non-monotonic behavior and changes 
sign as a function of K, the mean-field amplitudes are al- 
ways positive and monotonically increasing as function of 
k. As expected, the absence of fluctuations within MFT 
affects the quantitative estimate of the Casimir ampli- 
tude more strongly for the (o) BC than for the (+) BC. 



FIG. 28. Comparison of the scaling functions 6+(t, k) for 
a wall with a homogeneous (+) BC opposite to a chemically 
striped wall (Fig. [TJ as obtained for d = 3 and within MFT, 
i.e., for d = 4. The symbols are the data obtained from the 
MC simulations shown in Fig. 1131 The data obtained for 
k, — y oo [20J agree with the mean value of the data for (+, +) 
and (+, — ) BC of Ref. [29fl . The solid lines correspond to 
the MFT scaling functions #+ FT shown in Fig. [23]which have 
been rescaled according to Eq. (|64|l (see the main text and 
the caption of Fig. I13|) . Upon construction, for k — oo the 
positions and the heights of the maxima for d — 3 and d = 4 
agree. 




B. Scaling function of the critical Casimir force 

In order to compare also the temperature dependence 
of the scaling functions 9 + / q (t, k) of the critical Casimir 
force in d = 3 with their corresponding MFT estimates, 
it is useful to not only normalize the amplitude of the lat- 
ter but also to rescale them along the r axis by an over- 
all factor. Although this is an ad hoc procedure, it has 
turned out that a suitable combination of such rescaled 
MFT results with only partly available MC data might 
be a successful method in order to obtain quantitatively 
reliable approximations in an extended range of variables 
p7lj | . In the following we use a simple normalization of 
the MFT scaling functions 9^JJ{t,k). In Figs. [28] and 
I2U1 the mean-field scaling functions are rescaled linearly 



FIG. 29. Comparison of the scaling functions 9 (t, k) for a 
homogeneous (o) wall opposite to a striped wall (Fig. [2]). The 
symbols correspond to the MC data (d = 3) shown in Fig. 1201 
whereas the solid lines correspond to the MFT scaling func- 
tions (d = 4) shown in Fig. [24] which have been rescaled ac- 
cording to Eq. flM]). In contrast to Fig. [21] the rescaled MFT 
scaling functions differ qualitatively from the corresponding 
ones in d = 3. Upon construction, for k — oo the positions 
and the heights of the maxima for d = 3 and d = 4 agree. 
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so that for k -> oo the positions and the values of the 
maxima of the rescaled scaling functions O^JJ agree 
with those of the MC data. In Eq. (jMJ) T max ,+/o and 
7 mm+/o correspond to the position of the maximum of 
the scaling functions for k — > oo in d — 3 and d = 4, 
respectively. For the case of a homogeneous (+) wall 
opposite to a striped wall we can infer from the data 
of Ref. [29[ the rough estimates r maX!+ ~ —6.0 and 
#+(f m ax,+ , k — > oo) ~ 3.21 in d = 3 (see the caption 
of Fig. [Tl and Refs. SHI) and r^ T + ~ -31.960 and 
0^ FT (r^ + ,K -> oo) ~ 2.7531|6 (+ ; ) | in d = 4 (by 
taking the mean value of the scaling functions for (+, +) 
and (+,-) BC from Ref. [H|; see Fig. [23]). For a ho- 
mogeneous (o) wall opposite to a striped wall one has 
T max ,o = -1.174(10) and 6> (r maX!0 , k -> oo) = 0.564(3) 
in d = 3 (see Ref. [31| which agrees with the result shown 
in Fig.0 and ^ -7.0275 and 0^ FT (r^^ 

oo) ~ 0.35280|6 (+ . +) | in d = 4 as obtained from Ref. |2l| . 

Figure [28] shows the comparison of the scaling func- 
tions of the critical Casimir force for a homogeneous (+) 
wall opposite to a striped wall (see Fig. [1]). All MFT 
curves have been rescaled by the same factors according 
to Eq. ([64]) so that the position and the height of the 
maximum of the MFT curve for k — > oo agrees with the 
one obtained from the MC simulations in d = 3. As can 
be inferred from Fig. [2H1 the rescaled MFT behaviors as 
a function of r show a qualitative agreement with the 
corresponding MC results even for finite values of k. 

In Fig.[5n]we compare the scaling functions of the crit- 
ical Casimir force for a homogeneous wall with (o) BC 
opposite to a striped one (see Fig. [2]). The MFT scal- 
ing functions have been rescaled according to Eq. ([64"]). 
In contrast to the case shown in Fig. [25] these rescaled 
MFT scaling functions for the (o) case shown in Fig. [2j|] 
differ qualitatively from the corresponding behavior in 
d = 3. Whereas for k < 2 the MFT results suggest that 
the minima of the scaling functions exhibit a cusp-like 
singularity or a finite jump, the scaling functions 9 (t, k) 
in d = 3 are analytic at their minima. These differences 
are analogous to the ones obtained for homogeneous (o, o) 
BC at both surfaces [67| |68|j . 



VIII. SUMMARY AND OUTLOOK 

Within the Ising universality class we have studied 
the critical Casimir force for a film of thickness L by 
using Monte Carlo (MC) simulations in d — 3 spatial 
dimensions and by using mean-field theory. Along the 
lateral directions we have employed periodic boundary 
conditions (BC), whereas along the normal direction at 
the two confining surfaces fixed BC have been imposed. 
We have considered two cases: a homogeneous wall with 
(+) BC opposite to a wall patterned with alternating 
chemical stripes of equal width S + = 5_ with (+)/(—) 
BC (Fig. [1]) and a homogeneous wall corresponding to 
(o) BC opposite to a striped wall (Fig. [2]). In the limit of 



very narrow stripes, i.e., k = S+/L — > 0, the striped wall 
effectively mimics the behavior of Dirichlet (o) BC, so 
that for Kj — y the system reduces to the homogeneous 
cases with (+,o) or (o, o) BC, respectively (see Figs. [3] 
and 2]). In the opposite limit k — > oo, i.e., very broad 
stripes, in the first case (+; Fig. [1} the critical Casimir 
force equals the mean value of the corresponding forces 
for films with homogeneous (+,+) and (+,— ) boundary 
conditions at both surfaces, respectively. On the other 
hand, in the second case (o; Fig. [2]), deep in the two-phase 
regime, the corresponding limit is singular. 

We have investigated this system by combining MC 
simulations and numerical integration as well as by carry- 
ing out numerically the corresponding MFT calculation. 
We have employed an improved lattice model, for which 
the leading scaling corrections are suppressed. We have 
obtained the following main results: 

• In the finite size scaling limit the critical Casimir 
force Fc = L~ d d(T, k) per area and in units of fcgT 
is described (Eq. ©) by a universal scaling function 
8(t,k), with the scaling variables r = t{L/^) 1 / u 
and k = S+/L. Here t = {T-T c )/T c is the reduced 
temperature, £q is the non-universal amplitude of 
the correlation length £(i — > + ) = and 
S+ is the width of the stripes on the lower surface. 
In the limit k — > the patterned surface attains an 
effective Dirichlet BC (Eq. jl2j)). Within the range 
of aspect ratios p = L/L» (Figs. [THU) considered 
here, the MC data do not display a detectable de- 
pendence on p. Therefore we regard our results as 
the ones corresponding to the extrapolation to the 
film limit p — > 0. 

• In the limit of broad stripes, i.e., k ^> 1, the effects 
of the chemical steps separating the stripes vanish 
as oc k" 1 (Eqs. (O and ([11])). Thus, the total crit- 
ical Casimir force effectively approaches the sum of 
the forces between the individual stripes and the 
opposing wall. Accordingly, the assumption of ad- 
ditivity of the forces (which underlies the Derjaguin 
or proximity force approximation) generally holds 
for k — > oo. However, in the case of a homoge- 
neous wall with (o) BC opposite to a chemically 
striped wall, for k > 2 and r <C — 1, due to the 
formation of interfaces perpendicular to the film 
surfaces, the scaling function of the force varies as 
oc k -1 ^!^ oc L d /S + (for a fixed temperature t < 0; 
Eq. ([T5]) ). so that Fc does not decay for L — » oo 
as long as L < S+/2. Accordingly, for r — > — oo, in 
the subsequent limit k — > oo force additivity breaks 
down. The two limits k — > oo and r — > — oo do not 
commute. 

• By using MC simulations for d = 3, we have deter- 
mined the critical Casimir amplitude at T c for var- 
ious values of k, in the case of the BC illustrated in 
Figs. [T] and [2] as well as in the limit k — > 0, which 
corresponds to the BC shown in Figs. [3] and H] The 
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results are reported in Eqs. (|3~4"1) - ([3"T)) for the case of 
Fig. EQ and in Eqs. (|30 ]) -([25 f for the case of Fig. [2] 
Whereas in the first case involving a homogeneous 
(+) wall, the critical Casimir force is always repul- 
sive (Fig. [5]), in the case of a homogeneous (o) wall 
the critical Casimir amplitude is nonmonotonic and 
changes sign as a function of n (Fig. [6]). 

• Concerning T T c we have determined the criti- 
cal Casimir scaling functions & + / (t, k) in d = 3 for 
various values of well as in the limit k — > 0. In 
Figs.[7HT2land [TlHT51 we show the scaling functions 
8 + (t,k) and o (t, k), respectively, as determined 
for various film thicknesses. In Fig. [T3] we com- 
pare the universal scaling function # + (r, n) of the 
critical Casimir force between a homogeneous wall 
with (+) BC and a striped wall (Fig. [T|) for various 
values of determined from systems with the 
largest film thickness considered here, i.e., L z = 24, 
where L z = L/a + 1 and a is the MC lattice con- 
stant. We also compare our results with the uni- 
versal scaling function for the geometry consisting 
of a single chemical step (in the limit of vanishing 
aspect ratio studied in Ref. 20] ) which corresponds 
to the limit n — > oo. In this case, within the en- 
tire range < K < oo the critical Casimir force is 
always repulsive. 

• In contrast, for the case of a homogeneous (o) wall 
opposite to a striped one (Fig. [5]), the scaling func- 
tion of the critical Casimir force exhibits a rather 
different behavior. As shown in Fig. [501 the critical 
Casimir force varies nonmonotonically and changes 
sign as a function of k as well as a function of r. 
Moreover, for r < and for finite values of k the 
force may become much stronger than the ones for 
its limiting homogeneous cases (o, o) and (+, o) at- 
tained for k — y and k — > oo, respectively, which 
are also shown in Fig. [301 At k = 2 the system 
exhibits a transition of ground states from homo- 
geneous states for k < 2 to vertically striped states 
for k > 2. Whereas the scaling functions of the 
critical Casimir force for the homogeneous states 
exhibit a minimum at finite r < and vanish for 
t — > — oo, for k > 2 the scaling functions diverge 
for r — > — oo as oc \t\^ in accordance with Eq. (fT5j). 
This is confirmed by the MC results for k — 3 as 
shown in Fig. [TSJ Thus, the scaling functions of 
the critical Casimir force as obtained for this case 
- belonging to the Ising bulk universality class - 
do not vanish for r — > — oo. So far this peculiar 
feature is only known for the critical Casimir force 
acting in films belonging to the XY bulk univer- 
sality class and thus exhibiting Goldstone modes 

• In Scc. lVIl within MFT we have calculated the cor- 
responding scaling functions for the critical Casimir 
force for the two cases sketched in Fig.[T]and Fig. [3] 



The results for the suitably reduced critical Casimir 
amplitudes are shown in Figs. |2"T1 and l2"2l as a func- 
tion of n within a wide range of values. For k> 1 
the numerical MFT results agree with the asymp- 
totic behaviors of the scaling functions 0+ and Q a 
given in Eq. (fT3"]) and Eq. (fT4"|) . respectively, ac- 
cording to which they approach their correspond- 
ing limits for k — > oo as k . The suitably reduced 
universal scaling functions for r 7^ 0, as obtained 
within MFT, are shown in Fig.[33]for various values 
of k in the case of a homogeneous surface with (+) 
BC opposite to a striped surface. They interpolate 
smoothly between their limiting cases and always 
correspond to a repulsive critical Casimir force. 

• In the case of a homogeneous surface with (o) BC 
opposite to a striped surface the reduced MFT scal- 
ing functions are presented in Fig. 1241 They show 
a rich dependence on k. For k < 2 and r < 
the scaling function exhibits a minimum, and our 
numerical data suggest a cusp-like singularity or a 
finite jump of the scaling function at its minimum. 
For k > 2 the scaling functions diverge for r — > — 00 
and the MFT scaling functions agree to large ex- 
tent with the interface estimate given by Eq. ([15)) 
(Fig. US]). 

• The comparison of the suitably normalized Casimir 
amplitudes as obtained from MC simulations in 
d = 3 with the corresponding MFT ones reveals a 
good agreement for a homogeneous (+) surface op- 
posite to a striped surface (Fig. l2"o]l but qualitative 
differences for the corresponding (o) case (Fig. |2"T1) . 
Whereas in the latter case the data for d = 3 show 
a nonmonotonic behavior and a change of sign as 
function of k, in d = 4 the MFT Casimir ampli- 
tudes are always positive. 

• Similarly, as shown in Fig. [28j the behaviors of the 
full scaling functions 9+/ {t,k) as obtained from 
simulations and within MFT plus a suitable rescal- 
ing (Eq. (|64p ) agree qualitatively to large extent 
for a homogeneous (+) wall opposite to a striped 
one. On the other hand, for a homogeneous (o) 
surface opposite to a striped surface the MFT scal- 
ing functions show, even after rescaling, qualitative 
differences to the ones obtained via MC simulations 
(Fig.gSJ. However, within MFT as well as in d = 3, 
in the latter (o) case (Fig. [5]) we always observe a 
change of sign of the critical Casimir force from 
negative values at r <C — 1 to positive values for 
t > 0. In d = 3 this occurs for k > 3 and within 
MFT for all values of k > 0.5. At a fixed reduced 
temperature t, this zero at r = tq corresponds to 
a stable distance Lo(t) = Co"( r o/i) ! ' at which the 
upper plate levitates above the lower plate due to 
critical Casimir forces alone. The levitation height 
Loft) varies very sensitively as function of the re- 
duced temperature t = (T — T c )/T c . 
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• The computation of the critical Casimir force re- 
quires to subtract the bulk free energy density from 
the free energy density of the film. This bulk quan- 
tity is independent of the BC. For d — 3 we have 
determined it using a combination of MC simula- 
tions and numerical integration (see Add. |B|) . 

The present study is relevant for the critical behavior 
of films belonging to the Ising universality class and in 
the presence of a chemically structured substrate. This 
can be experimentally realized by considering complete 
wetting films of binary liquid mixtures near their criti- 
cal endpoints of demixing and by exposing their vapor 
phases to a chemically structured substrate [l(| EH ■ The 
critical Casimir forces can be inferred by monitoring the 
thicknesses of the wetting films. This realizes the (+) BC 
versus a striped surface. The surface fields describe the 
preferences of the two species for the confining interfaces 
of the wetting films. 

Another realization consists of studying directly the 
force acting on a colloidal particle immersed in a critical 
binary liquid mixture and exposed to a chemically struc- 
tured substrate, as has been done in Ref. [T3|- In this 
case the normal critical Casimir force is approximately 
the one for the film geometry investigated here, provided 
the radius of the colloidal particle is sufficiently large rel- 
ative to its distance from the wall. However, near T c 
for such a system an additional lateral critical Casimir 
force sets in. In Ref. [4l| the critical Casimir force for a 
sphere in front of a chemically structured substrate has 
been studied by means of mean-field theory as well as 
in d = 3 by using the Derjaguin approximation. In this 
study, it was found that for suitable geometric features of 
the stripes on the substrate and in the presence of homo- 
geneous BC on the spherical colloid levitation is possible 
even for r > 0, i.e., in the homogeneous phase of the 
solvent. Although these experimental studies [THHEI are 
closely related to the setup studied here, a re-evaluation 
of the existing data is not sufficient in order to compare 
them with the present theoretical predictions. On one 
hand, the authors of Refs. [3, [l5[ have measured only 
the lateral forces acting on the colloidal particles, and 
not the normal ones studied here. On the other hand, 
in order to effectively mimic the film geometry studied 
here, the radius of the colloidal particles should be much 
larger than the stripe widths 5+ and S-, whereas the 
length scales realized in those experimental studies are 
of the same size [2 UU . 

In view of recent MC results for the critical Casimir 
force of a sphere in front of a homogeneous wall [34j, it 
would be very interesting to extend this study by con- 
sidering a sphere in front of a chemically structured wall. 
Besides analyzing directly the latter experimental setup, 
this would also provide the possibility of elucidating the 
range of validity of the Derjaguin approximation, which 
is commonly employed for curved geometries (40l . l4ll ] . 

Here we have determined the critical Casimir force in 
the presence of a chemically structured substrate by us- 
ing MC simulations in spatial dimension d = 3, and by 



using mean-field theory, which holds in d — 4. In order to 
complement this spatial dependence and to further probe 
the relevance of fluctuations, it would be interesting to 
investigate the corresponding system in d = 2, where 
some exact results are available [43j and conformal in- 
variance allows one to determine exactly certain critical 
properties. 

The present study also lends itself to further exten- 
sions. Here we have considered stripes with (+) and (— ) 
BC of equal widths. A natural extension of the present 
study would consist of calculating the critical Casimir 
force as a function of the ratio of the widths of the (+) 
and (— ) stripes. Moreover, by considering two striped 
surfaces, one can also investigate the corresponding lat- 
eral critical Casimir force. So far the case of two striped 
surfaces has been investigated by mean-field theory for 
the film geometry [38| ; the issue of the lateral force has 
been analyzed by mean-field theory and in d = 3 within 
the Derjaguin approximation for the sphere- wall geome- 
try [HI, 

Finally, as mentioned in Sees. IVTlVIl for the BC shown 
in Fig. [5] and for k — 2 the system displays a rich glassy 
behavior at low temperatures. This deserves further in- 
vestigations. 
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Appendix A: Monte Carlo simulations 

In this appendix we report certain technical details of 
the MC simulations we have performed. As explained in 
Sec. lIVI the evaluation of the Casimir force at T c has been 
carried out in two steps. First, we have determined the 
thermal average (H2 — %i)a which appears in Eq. ([2"5]l, 
This is obtained by a standard MC simulation for the 
ensemble characterized by the crossover Hamiltonian H\ 
defined in Eq. (|25|) . We have implemented a combination 
of the standard Metropolis and Wolff cluster algorithms: 
each MC step consists of 1 Metropolis sweep over the 
entire lattice in lexicographic order and L z Wolff single- 
cluster flips; L z denotes the total number of lattice layers 
including the surfaces of fixed spins, so that there are 
L z — 2 layers of fluctuating spins in the case of the BC 
shown in Fig. [1] L z — 1 layers in the case of BC shown 
in Figs. [2] and |3l and L z layers in the case of the (o, o) 
BC illustrated in Fig. 2) As random number generator 
we have used the double precision SIMD-oriented Fast 
Mersenne Twister (dSFMT) [zl]. Important details of 
the simulations performed at the critical temperature are 
reported in Tables fVTlIXI Additional details concerning 
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the implementation of the simulating algorithm can be 
found in Ref. [!}. 

As explained in Sec. El the determination of the 
scaling function for the critical Casimir force has been 
obtained by sampling the reduced energy densities 
E(p', L z ,L X7 s+) and E(p', L z - 1, L x , a+) (see Eq. (JMH) 
followed by carrying out numerically the integration in 
Eq. (IM)) by using Simpson's rule. An upper bound of the 
systematic error due to the discretization of the integrals 
can be determined by sampling the fourth derivative of 
the integrand: by computing d 4 E((3, L z , L x , s + )/ (df3 4 ) 
we have checked that such a systematic error is always 
negligible compared to the statistical errors. (Since for 
L z — >• oo the quantity d 4 E(f3 7 L z , L x , s + )/ (df3 4 ) diverges 
at the critical point, the number of sampled points has to 
increase with L z .) In TablclXlwe report important details 
concerning these simulations associated with Eq. (f5~4"|). 
For each film thickness and BC we have considered the 
same three aspect ratios p for determining the scaling 
functions as the ones used for determining the critical 
Casimir amplitude (see Tables IVllIX[) . For the BC shown 
in Figs. [T] and O we have verified that the sampled re- 
duced energy densities are de facto independent of p. 
Therefore our results capture reliably the limit p — > 0; 
we have averaged them over the three aspect ratios con- 
sidered. Concerning the BC shown in Figs. [2] and HI as 
discussed in Sec. [V] the data exhibit a weak dependence 
on the aspect ratio at low temperatures and we have con- 
sidered the three aspect ratios separately, i.e., without 
taking this average. 

Finally, we mention that with the above described sim- 
ulation algorithm and for the BC shown in Fig[3J we oc- 
casionally observed the appearance of metastable states 
at low temperatures, which cause the thermalization of 
the run to be rather long. We have found that this prob- 
lem can be healed by starting the simulations with an 
ordered states. 



Appendix B: Determination of the bulk free energy 
density 

Here we report certain details concerning the determi- 
nation of the bulk free energy density which is needed for 
calculating the critical Casimir force (see Eqs. (|50p and 
(|53[) ) . For this purpose we have simulated the improved 
Blume-Capel model described by Eq. (fTB")) for a simple cu- 
bic lattice with periodic BC in all directions and lattice 
sizes L z — 24 — 256. For this system we have determined 
the reduced energy density E(/3, L z ) and the reduced free 
energy density F(j3, L z ) as defined in Eqs. (|2"tj|) and (|2"T1) . 
For the sake of simplicity, here we omit the dependence 
on L x and s + because the lattice considered here has the 
same size in all directions and it does not have any sur- 
face. Since the aim is to extract the thermodynamic limit 
of these quantities from finite-size results, we recall the 
expected behavior of the corresponding finite-size parts. 
For T y^T c and L z ^> £, E(f3, L z ) approaches its infinite- 



volume limit -Ebuik(/3) as 

6E(0, L z ) = E{0, L z ) - £ b ulk(/3) ~ {L z /£) k+l e~ L ^, 

(Bl) 

where k is an integer. Conversely, in the region where 
£ « L z , one has (a — 2 — iv) 



SE(f3,L z )=t 1 - a h E (L,/f) = 



1 



h E (L z /0, (B2) 



where the scaling function 1ie(x) is universal up to a 
prefactor and /is(x) = O(l) for £ s» L z . The reduced 
free energy density F(j3, L z ) can be obtained by integrat- 
ing E((3,L Z ) according to Eq. (|2"2"j) . It follows that, for 
T > T c and L z ^> £, F(f3,L z ) approaches its infinite- 
volume limit -Fbuik(/3) as 



SF(J3,L, 



F{p, L z ) - F bulk (/3) ~ (L z /0 k e- L =/Z. 



(B3) 

In deriving Eq. (|B3I) . we have used the fact that for 
T > T c the condition L z 3> £ is satisfied throughout 
the interval of integration on the r.h.s. of Eq. ((22)) . This 
is not the case if T < T c . For T < T c , L z > £, and by 
using Eq. (|22l) . the finite-size correction 6F(j3, L z ) can be 
expressed as 

5F(0,L z ) = F{p,L z )-F hulk (0) 

rP (B4) 
= 5F{P -> oo, L z ) + / dp'5E{P',L z ). 

J oo 

In the second term of the r.h.s. of Eq. (IB4|) one has 
L z S> £ throughout the integration interval. Thus, by us- 
ing Eq. (|B1[) . the integral on the r.h.s. of Eq. (|B4I) varies 
as (L z /^) k e~ Lz ^. The finite-size correction SF((3 — > 
oo,L z ) can be inferred from computing F(/3,L Z ) for 
(3 — s- oo and for a finite size L z . For j3 — > oo, the Gibbs 
measure is dominated by the twofold degenerate ground 
state, consisting of a configuration in which all spins are 
fixed to +1 or to —1. By using the definition of F(f3,L z ) 
given in Eq. ([2"Tj) , one has 



F(/3^oo,L 2 
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(B5) 



where D is the coupling constant appearing in the second 
term of the Hamiltonian given in Eq. (flT)}) . By taking the 
L z — > oo in Eq. (|B5I) . we identify the second term on the 
r.h.s. of Eq. (|B5|) as the infinite-volume limit i 7 buik(/3) 



7g. Thus we infer 5F{fi -> oo,L z ) = (\n2)/(L 3 z ). Thus, 
for T <T C and L z ^> £, F(f3, L z ) approaches its infinite- 
volume limit fbuikCS) as 

In 2 



6F(P,L z )~(L z /t) k e- L ^ + 



(B6) 



From Eqs. (J22J) and (JB6J) one finds that 



6E(/3',L z )d/3' 
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TABLE V. The total number N a t eps of MC steps and the number iVth C rm of MC steps discarded in order to achieve thermalization 
as used to determine the critical Casimir amplitudes for film thicknesses L z > 24, for aspect ratios p = L z /L x < 1/8, and for 
the BC shown in Figs. [T] and Every MC step consists of 1 Metropolis sweep over the entire lattice and L z Wolff single-cluster 
flips. Additional details concerning the simulation algorithm can be found in Ref. [20|. 
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350 


70 


48 


1/16 


350 


50 


48 


1/16 


350 


50 


48 


1/24 


170 


34 


TABLE VI. Same as Table|V]for k = 1, 2, 3. 






k — y : 


(o,o) 






K — 


1/4 






ft = 1/2 




L z 


P 


AWs/10 3 


iVtherm/10 3 


L z 


P 


Astops/lO 1 


3 AT therm /10 3 


L z 


P 


iVsteps/10 3 


JVtherm/10 3 


24 


1/8 


12000 


1200 


24 


1/8 


12000 


1200 


24 


1/8 


12000 


1200 


24 


1/12 


8000 


800 


24 


1/12 


8000 


800 


24 


1/12 


8000 


800 


24 


1/16 


6000 


600 


24 


1/16 


6000 


600 


24 


1/16 


6000 


600 


32 


1/8 


32000 


600 


32 


1/8 


6000 


600 


32 


1/8 


6000 


600 


32 


1/12 


16000 


300 


32 


1/12 


3000 


300 


32 


1/12 


3000 


300 


32 


1/16 


8000 


150 


32 


1/16 


1500 


150 


32 


1/16 


1500 


150 



TABLE VII. Same as Table |V] for the BC of Figs. Hand g] 



ft = 3/4 ft = l 



L z 


P 


iVstcps/10 3 


iV therm /10 3 


L z 


P 


iVstcps/10 3 


AWrm/10 3 


24 


1/9 


12000 


120 


24 


1/8 


12000 


1200 


24 


1/12 


8000 


800 


24 


1/12 


8000 


800 


24 


1/15 


6000 


600 


24 


1/16 


6000 


600 


32 


1/9 


6000 


600 


32 


1/8 


6000 


600 


32 


1/12 


3000 


300 


32 


1/12 


3000 


300 


32 


1/15 


1500 


150 


32 


1/16 


1500 


150 



TABLE VIII. Same as Table [VII] for ft = 3/4 and 1. 
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K — 


2 






k — 3 






P 


AUps/10 3 


AWrm/10 3 




P 


JVsteps/10 3 


Wtherm/lO 3 


24 


1/24 


600 


100 


24 


1/24 


600 


60 


24 


1/36 


400 


80 


24 


1/36 


400 


60 


24 


1/48 


300 


60 


24 


1/48 


300 


60 


32 


1/24 


150 


30 


32 


1/24 


150 


20 


32 


1/36 


70 


15 


32 


1/36 


70 


10 


32 


1/48 


40 


8 


32 


1/48 


40 


8 



TABLE IX. Same as Table [VTT] for k = 2 and 3. 



L z p p max A/3 



8 


0.327721735 


0.427721735 


0.0005 


12 


0.327721735 


0.427721735 


0.0001 


16 


0.327721735 


0.427721735 


0.0005 


24 


0.377721735 


0.397721735 


0.00002 



TABLE X. The lowest (/3 ) and the highest (/3 max ) inverse 
temperatures used for the computation of the scaling func- 
tions associated with the free energy differences via Eq. (|54[) . 
The integrals have been computed numerically using Simp- 
son's rule, with the reported intervals A/3 between two con- 
secutive points. For each film thickness we have considered 
the same three aspect ratios as the ones used for determining 
the critical Casimir amplitude (see Tables [VllXfl . 

where the support of the integrand is actually confined 
to the region where L z w £. In this region the finite size 
correction of the reduced free energy density is given by 

8F(p,L M ) = jph F (L a /t), (B8) 

where, as in Eq. (IB2[) . the universal scaling function 
}if(x) = O(l) for £ Ri L z . A comparison of the finite- 
size corrections for the reduced energy density given in 
Eqs. (|Bll) - (|B2p and those for the reduced free energy den- 
sity in Eqs. (|H3]l . ([B"6j). and §gf\ shows that F(P,L Z ) 
converges faster to limit for L z — > oo than E(f3, L z ). The 
only exception to this rule occurs in the low-temperature 
phase, T < T c and L z >• £, where the reduced free en- 
ergy density exhibits an additional finite-size correction 
In2/L 3 (see Eq. (|B7|) L However, because this correction 
term is known exactly, one can eliminate it by subtract- 
ing it explicitly. 

In order to compute the bulk free energy density, we 
proceed as follows. At a given lattice size L z , we com- 
pute the reduced energy density E(j5,L z ) in an interval 
[Anin, Anax] around the inverse critical temperature j5 c = 
0.387721735(25) [45|. In order to minimize the error bars 
we have implemented the control- variates scheme intro- 
duced in Ref . [75[ . Control variates are observables which 
have a vanishing mean value and therefore can be added 
to any observable without changing its mean value; con- 
trol variates provide also an additional check of the MC 
simulations. In the second step F(P,L Z ) — F(/3 min , L z ) 



L 


Pmi n 


/^m ax 


A/3 


24 


0.327721735 


0.427721735 


0.0002 


32 


0.347721735 


0.427721735 


0.0001 


48 


0.367721735 


0.407721735 


0.0001 


64 


0.377721735 


0.397721735 


0.0001 


96 


0.380521735 


0.395721735 


0.00005 


128 


0.381521735 


0.394521735 


0.00005 


192 


0.384121735 


0.393321735 


0.00002 


256 


0.385321735 


0.391321735 


0.00001 



TABLE XI. The interval of integration [/3 m i n , /3 max ] for each 
lattice size L used in the determination of the bulk free en- 
ergy density. We have implemented Simpson's rule with the 
reported distances A/3 between two consecutive points. 



is calculated by numerically integrating Eq. (|22p . For 
this purpose we have used Simpson's rule. The resulting 
quantity F(f3,L z ) — F(/3 m j n , Lz) suffers from two types 
of errors: a statistical error originating from the statisti- 
cal error bars of the integrand E((3, L z ) and a systematic 
error due to the chosen quadrature. In the present case 
and as mentioned above, the maximum systematic er- 
ror in Simpson's rule can be computed by estimating the 
fourth derivative of E(P,L Z ). We have always checked 
that such an error is at least one order of magnitude 
smaller than the statistical error, so that it can be safely 
neglected and the statistical error bar is a correct mea- 
sure of the uncertainty of the reduced free energy den- 
sity. The integration of E(j3,L z ) leads to the value of 
F{P,L Z ) — F(f3 m i D , L z ) for several inverse temperatures 
P € [/3 m in, Anax]- For those values of P < fi c for which 
L z ^> £, we regard our results for finite L z to be the ones 
for infinite L z if the statistical error bars are smaller than 
the finite-size correction. To this end, we have checked 
that E(P,L Z ) is, within the numerical accuracy, inde- 
pendent of L z by comparing the values obtained for two 
consecutive lattice sizes: as discussed above, E(P,L Z ) is 
expected to converge to the thermodynamic limit slower 
than F(P,L Z ). Roughly speaking, with the present nu- 
merical accuracy, the finite-size scaling corrections are 
negligible for L z /t; < 20. For T < T c we use the more 
conservative bound L z /t; < 35—40, and we explicitly sub- 
tract the additional finite-size term (In2)/L 3 which ap- 
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pears in Eq. (IB7|) . We note that the non- universal ampli- 
tude £Z of the correlation length below T c is roughly half 
°f £02 : £oi/£o2 = 1-957(7) (76[. At any given lattice size 
L z , this procedure results in the estimate of the bulk free 
energy density for a subset [/3 min , /3 inf ] U [f3 sup , /3 max ] of the 
integration interval [/3 min , /3 max ], with £(Anf) ~ L z /20 
and ^(Aup) « L,/40. Thus, for /3 < /3 lnf and /3 > /3 8up , 
F((3,L Z ) — F((3 min ,L z ) agrees within error bars with 
-Fbuik(/3) - fbuikCAmn), while [/3 inf , /3 sup ] is the interval in 
which the finite-size correction 5F(/3,L Z ) is not negligi- 
ble. In the next step we have applied the above procedure 
for a larger lattice size L' z > L z and the smaller integra- 
tion interval [/3' min = /3 ini , j3' max = /3 sup ]. This results 
in the quantity F(J3, L' z ) — F(f3' min , L' z ) to which we add 

^XAnin' Lz) — F((3 min , L z ) ~ i 7 'bulk(/3min) — ^bulk(/3min) as 

determined from the lattice size L z , so that we finally 
obtain the desired quantity F((3,L' Z ) — i 7 buik(/3min)- As 
before, this results in the estimate of the bulk free en- 
ergy for /? G [/3 mi „, (3' in{ ] U [^ up ,/3 max ], with # nf > ft nf 
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